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Abstract 

Assuming that the observed pattern of 3-neutrino mixing is related to the existence of a 
(lepton) flavour symmetry, corresponding to a non-Abelian discrete symmetry group Gf, 
and that Gf is broken to specific residual symmetries G e and G„ of the charged lepton 
and neutrino mass terms, we derive sum rules for the cosine of the Dirac phase <5 of the 
neutrino mixing matrix U. The residual symmetries considered are: i) G e = Z 2 and 
G v = Z n , n > 2 or Z n x Z m , n,m > 2; ii) G e = Z n , n > 2 or Z n x Z rn , n,m > 2 and 
G u = Z 2 ; iii) G e = Z 2 and G u = Z 2 , iv) G e is fully broken and G„ = Z n , n > 2 or 
Z n x Z m , n, m > 2; and v) G e = Z n , n > 2 or Z n x Z m , n, m > 2 and G v is fully broken. 
For given G e and G v , the sum rules for cos <5 thus derived are exact, within the approach 
employed, and are valid, in particular, for any Gf containing G e and G„ as subgroups. 
We identify the cases when the value of cos <5 cannot be determined, or cannot be uniquely 
determined, without making additional assumptions on unconstrained parameters. In a 
large class of cases considered the value of cos 6 can be unambiguously predicted once 
the flavour symmetry G/ is fixed. We present predictions for cos 5 in these cases for the 
flavour symmetry groups G/ = S 4 , A 4 , T' and A$, requiring that the measured values of 
the 3-neutrino mixing parameters sin 1 2 6 * 12 , sin 2 613 and sin 2 023, taking into account their 
respective 3 a uncertainties, are successfully reproduced. 
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1 Introduction 


The discrete symmetry approach to understanding the observed pattern of 3-neutrino mixing 
(see, e.g., 0 ) , which is widely explored at present (see, e.g., 00 ) , leads to specific correlations 
between the values of at least some of the mixing angles of the Pontecorvo, Maki, Nakagawa, 
Sakata (PMNS) neutrino mixing matrix U and, either to specific fixed trivial or maximal 
values of the CP violation (CPV) phases present in U (see, e.g., (&-10| and references quoted 
therein), or to a correlation between the values of the neutrino mixing angles and of the Dirac 
CPV phase of U 11 15 d] As a consequence of this correlation the cosine of the Dirac CPV 


phase 5 of the PMNS matrix U can be expressed in terms of the three neutrino mixing angles 
of U 11-14 , i.e., one obtains a sum rule for cos 5 . This sum rule depends on the underlying 


discrete symmetry used to derive the observed pattern of neutrino mixing and on the type of 
breaking of the symmetry necessary to reproduce the measured values of the neutrino mixing 
angles. It depends also on the assumed status of the CP symmetry before the breaking of the 
underlying discrete symmetry. 

The approach of interest is based on the assumption of the existence at some energy scale 
of a (lepton) flavour symmetry corresponding to a non-Abelian discrete group Gf. Groups 
that have been considered in the literature include S4, A4, T ', A 5 , D n (with n = 10,12) 
and A( 6 n 1 2 ), to name several. The choice of these groups is related to the fact that they 
lead to values of the neutrino mixing angles, which can differ from the measured values at 
most by subleading perturbative corrections. For instance, the groups A4, S4 and T' are 
commonly utilised to generate tri-bimaximal (TBM) mixing 18 ; the group S4 can also be 
used to generate bimaximal (BM) mixing 19 10 A 5 can be utilised to generate golden ratio 
type A (GRA) 121—231 mixing; and the groups Dio and D \2 can lead to golden ratio type B 
and hexagonal (HG) 125 mixing. 


(GRB) 124] 

The flavour symmetry group Gf can be broken, in general, to different symmetry sub¬ 
groups G e and G u of the charged lepton and neutrino mass terms, respectively. G e and G v 
are usually called “residual symmetries” of the charged lepton and neutrino mass matrices. 
Given Gf , which is usually assumed to be discrete, typically there are more than one (but still 
a finite number of) possible residual symmetries G e and G u . The subgroup G e , in particular, 
can be trivial, i.e., Gf can be completely broken in the process of generation of the charged 
lepton mass term. 

The residual symmetries can constrain the forms of the 3x3 unitary matrices U e and U u , 
which diagonalise the charged lepton and neutrino mass matrices, and the product of which 
represents the PMNS matrix: 

u = u\u v . ( 1 ) 


Thus, by constraining the form of the matrices U e and U v , the residual symmetries constrain 
also the form of the PMNS matrix U. 

In general, there are two cases of residual symmetry G v for the neutrino Majorana mass 
term when a portion of Gf is left unbroken in the neutrino sector. They characterise two 
possible approaches — direct and semi-direct |4| — in making predictions for the neutrino 
mixing observables using discrete flavour symmetries: G u can either be a Z2 x Z2 symmetry 


1 In the case of massive neutrinos being Majorana particles one can obtain under specific conditions also 
correlations between the values of the two Majorana CPV phases present in the neutrino mixing matrix [16] 
and of the three neutrino mixing angles and of the Dirac CPV phase 11.17 . 

2 Bimaximal mixing can also be a consequence of the conservation of the lepton charge Lf = L e — — L T 

(LC) 


20 


supplemented by a H — r symmetry. 
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(which sometimes is identified in the literature with the Klein four group), or a Z 2 symmetry. 
In models based on the semi-direct approach, where G u = Z 2 , the matrix U u contains two free 
parameters, i.e., one angle and one phase, as long as the neutrino Majorana mass term does 
not have additional “accidental” symmetries, e.g., the fi — r symmetry. In such a case as well 
as in the case of G u = Z 2 x Z 2 , the matrix U v is completely determined by symmetries up to 
re-phasing on the right and permutations of columns. The latter can be fixed by considering 
a specific model. It is also important to note here that in this approach Majorana phases are 
undetermined. 

In the general case of absence of constraints, the PMNS matrix can be parametrised in 
terms of the parameters of U e and U v as follows |26|: 

U = UlU u = (U e )^U u Q 0 . (2) 

Here U e and U v are CKM-like 3x3 unitary matrices and \k and Qq are given by: 

m = diag ^1,e -1 ^, e - *^ , Q 0 = diag , (3) 


where i/j, ui , £21 and £31 are phases which contribute to physical CPV phases. Thus, in general, 
each of the two phase matrices and Q 0 contain two physical CPV phases. The phases in 
Qo contribute to the Majorana phases [l 6 | in the PMNS matrix (see further) and can appear 
in ecp ([ 2 ]) as a result of the diagonalisation of the neutrino Majorana mass term, while the 
phases in T can result from the charged lepton sector (u\ = (U e )^ ^), from the neutrino 
sector {U v = 'S>U u Qo), or can receive contributions from both sectors. 

As is well known, the 3x3 unitary PMNS matrix U can be parametrised in terms of 
three neutrino mixing angles and, depending on whether the massive neutrinos are Dirac or 
Majorana particles, by one Dirac CPV phase, or by one Dirac and two Majorana 
phases: 


16 CPV 


U = U\ U v = VQ , Q = diag (1, 


a 21 
2 , 


( 4 ) 


where 0:21,31 are the two Majorana CPV phases and V is a CKM-like matrix. In the standard 
parametrisation of the PMNS matrix [I], which we are going to use in what follows, V has 
the form: 



( 

C 12 C 13 

S12C13 

si 3 e- i5 

V = 


— S12C23 — Cl2'S23'Sl3e*' 5 

C12C23 - Si2S23Sise l5 

S 23 C 13 


V 

•S12S23 - Ci2C23Si3e l<5 

— C12S23 — Sl2C23Sl3e i<5 

C 23 C 13 


where 0 < 5 < 27t is the Dirac CPV phase and we have used the standard notation Cij = cos 9ij , 
Sij = sin Oij with 0 < 9^ < 7 r/ 2 . Notice that if CP invariance holds, then we have 5 = 0, n, 2ir, 
with the values 0 and 27 t being physically indistinguishable, and 021 = kir, 031 = k! tt, 
k,k' = 0, 1,2^] Therefore, the neutrino mixing observables are the three mixing angles, 612 , 
013, 023 1 the Dirac phase 5 and, if the massive neutrinos are Majorana particles, the Majorana 
phases 021 and 031 . 

The neutrino mixing parameters sin 2 0 i 2 , sin 2 023 and sin 2 0 i 3 , which will be relevant for 
our further discussion, have been determined with a relatively good precision in the recent 


3 In the case of the type I seesaw mechanism of neutrino mass generation the range in which 021 and 031 
Thus, in this case 021 and 031 possess CP-conserving values for k, k' = 0,1, 2, 3,4. 


vary is [0, 4 - 7 r] 


27 
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global analyses of the neutrino oscillation data 28,29 . For the best fit values and the 3er 
allowed ranges of sin 2 $i 2 , sin 2 623 and sin 2 # 13 , the authors of ref. 28 have obtained: 

(sin 2 6»i 2 )bf = 0.308 , 0.259 < sin 2 0 12 < 0.359 , 


(sin 2 (9 23 )bf = 0.437 (0.455), 0.374 (0.380) < sin 2 0 23 < 0.626 (0.641), 

(sin 2 0 i 3 ) bf = 0.0234 (0.0240), 0.0176 (0.0178) < sin 2 0 13 < 0.0295 (0.0298). 


( 6 ) 

( 7 ) 

( 8 ) 


Here the values (values in brackets) correspond to neutrino mass spectrum with normal or¬ 
dering (inverted ordering) (see, e.g., |1|), denoted further as the NO (IO) spectrum. 

In ref. |ll] (see also 12-14]) we have considered the cases when, as a consequence of 
underlying and residual symmetries, the matrix U u , and more specifically, the matrix U v in 
eq. ([2^, has the i) TBM, ii) BM, iii) GRA, iv) GRB and v) HG forms. For all these forms we 
have U y = i? 23 ( 02 3 ) 7 ?i 2 ( 0 i 2 ) with 0 23 = — tt/ 4, i ? 23 and R 12 being 3x3 orthogonal matrices 
describing rotations in the 2-3 and 1-2 planes: 


Uu = R-23 (0 U 2 3 ) 7 ? 12 (0 U 12 ) = 


/ 

COS 0 \2 

sin 0/2 

° ^ 


sin 0\ 2 

cos 0'( 2 

1 


V2 

V2 

V2 


sin d'{ 2 

COS 0y 2 

1 

V 

V2 

V2 

v/2 / 


( 9 ) 


The value of the angle 0 32 , and thus of sin 2 0\ 2 . depends on the form of U v . For the TBM, 
BM, GRA, GRB and HG forms we have: i) sin 2 0 / 2 = 1/3 (TBM), ii) sin 2 0 / 2 = 1/2 (BM), 
iii) sin 2 = (2 + r) -1 = 0.276 (GRA), r being the golden ratio, r = (1 + y/b)/2, iv) 
sin 2 0 \ 2 = (3 - r)/4 = 0.345 (GRB), and v) sin 2 0 \ 2 = 1/4 (HG). 

The TBM form of U u , for example, can be obtained from a Gf = A 4 symmetry, when 
the residual symmetry is G v = Z 2 . i.e. the S generator of A 4 (see Appendix |A|) is unbroken. 
In this case there is an additional accidental fi — r symmetry, which together with the Z 2 
symmetry leads to the TBM form of U v (see, e.g., { 3 ]). The TBM form can also be derived from 
Gf = T' with G y = Z 2 , provided that the left-handed (LH) charged leptons and neutrinos each 
transform as triplets of T' and the TST 2 element of T' is unbroken, see Appendix [A|for further 
explanation. Indeed when working with 3-dimensional and 1-dimensional representations of 
T' , there is no way to distinguish T' from A 4 


30 
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and listed above we have i) 0/ 3 = 0 , which should 


Finally, one can obtain BM mixing from, 
e.g., the Gf = S 4 symmetry, when the residual symmetry is G u = Z 2 . There is an accidental 
H — t symmetry in this case as well |31 
For all the forms of U u considered in 

be corrected to the measured value of 0i 3 = 0.15, and ii) sin 2 0 23 = 0.5, which might also need 
to be corrected if it is firmly established that sin 2 0 23 deviates significantly from 0.5. In the 
case of the BM and HG forms, the values of sin 2 0 \ 2 lie outside the current 3cr allowed ranges 
of sin 2 0\ 2 and have also to be corrected. 

The requisite corrections are provided by the matrix U e , or equivalently, by U e . The 
approach followed in 11-14 corresponds to the case of a trivial subgroup G e , i.e., of Gf com¬ 


pletely broken by the charged lepton mass term. In this case the matrix U e is unconstrained 


and was chosen in 11 on phenomenological grounds to have the following two forms: 

U e = R/ 2 \0t 2 ), 


Ue = R 2 i{0 e 23)R/2^l 2 
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( 10 ) 

( 11 ) 


















These two forms appear in a large class of theoretical models of flavour and theoretical studies, 
in which the generation of charged lepton masses is an integral part (see, e.g., 0 32-{37]). 

In this setting with U u having one of the five symmetry forms, TBM, BM, GRA, GRB 
and HG, and U e given by eq. |ll|), the Dirac phase 5 of the PMNS matrix was shown in 111 
to satisfy the following sum rulej_j 


cos 5 = 


tan 023 


sin 2 #i 2 sin 0 i 3 


[cos 2012 + (sin 2 0i2 — cos 2 0i 2 ) (l — cot 2 023 sin 2 013)] . (12) 


Within the approach employed this sum rule is exact 0 It is valid, in particular, for any value 
of the angle 0 23 14 In 


11 


by using the sum rule in eq. (12), predictions for cos <5 and S 
were obtained in the TBM, BM, GRA, GRB and HG cases for the best fit values of sin 2 0i2, 
sin 2 023 and sin 2 0i3. The results thus obtained permitted to conclude that a sufficiently 
precise measurement of cos 5 would allow to discriminate between the different forms of U u 
considered. 


Statistical analyses of predictions of the sum rule given in eq. (12) i) for 5 and for the Jcp 


factor, which determines the magnitude of CP-violating effects in neutrino oscillations 38 
using the current uncertainties in the determination of sin 2 0 i 2 , sin 2 0 i 3 , sin 2 023 and <5 from 


28 , and ii) for cos <5 using the prospective uncertainties on sin 2 0 i 2 , sin 2 0 i 3 and sin 2 023 , were 


of 


performed in [13| for the five symmetry forms — BM (LC), TBM, GRA, GRB and HG 
Uu- 

In |14j we extended the analyses performed in |ll|l3| by obtaining sum rules for cos 5 for 
the following forms of the matrices U e and Uu^\ 

A. U u = ^ 23 (^ 23 )^ 12 (^ 12 ) with $23 = —tt/ 4 and 0^ 2 as dictated by TBM, BM, GRA, 


GRB or HG mixing. 
Rn^h)Rn(0f 2 y, 


and i) Ue = R^(0 e 13 ), ii) Ue = R 23 \eh)Rii ($ 13 ), and hi) U e = 


B. U v = ^ 23 (^ 23 )^ 13 (^ 13 )^ 12 (^ 12 ) with 023) $13 an d 0i 2 fixed by arguments associated with 
symmetries, and iv) U e = R^^O 12 ): and v) U e = 

The sum rules for cos <5 were derived first for 0 23 = — 7r/4 for the cases listed in point A, and 
for the specific values of (some of) the angles in U v , characterising the cases listed in point 
B, as well as for arbitrary fixed values of all angles contained in U v . Predictions for cos <5 
and Jcp (cos<5) were also obtained performing statistical analyses utilising the current (the 
prospective) uncertainties in the determination of sin 2 0i2, sin 2 0i3, sin 2 023 and 5 (sin 2 0i2, 
sin 2 0 i 3 and sin 2 623 )- 

In the present article we extend the analyses performed in [IP 13,[l4j to the following 


cases: 


1) G e = Z2 and G v = Z n , n > 2 or Z n x Z m , n, m > 2; 

2) G e = Z n , n > 2 or Z n x Z m , n, m > 2 and G u = Z2; 


The sum rule is given in the standard parametrisation of the PMNS matrix (see, e.g., 1 ). 

For the TBM and BM forms of U v , and for U e given in eq. (Ill, it was first derived in ref. 12 . 


°The two forms of U e in eqs. (10 1 and (jllj) lead, in particular, to different predictions for sin 023: for 
023 = — 7 t/4 in the case of eq. (101 we have sin“ 023 = 0.5(1 — sin 2 013 ), while if U e is given by eq. (Ill, sin 2 023 
can deviate significantly from 0.5. 


We performed in 14 a systematic analysis of the forms of U e and Lb, for which sum rules for cos S of the 


type of eq. (121 could be derived, but did not exist in the literature. 
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3) G e = Z '2 and G u = Z 2 ; 

4) G e is fully broken and G u = Z n , n > 2 or Z n x Z m , n,m > 2; 

5) G e = Z n , n > 2 or x Z m , n,m > 2 and is fully broken. 


In the case of G e = Z 2 (G„ = Z 2 ) the matrix U e (U u ) is determined up to a U{2) transfor¬ 
mation in the degenerate subspace, since the representation matrix of the generator of the 
residual symmetry has degenerate eigenvalues. On the contrary, when the residual symmetry 
is large enough, namely, G e = Z n , n > 2 or Z n x Z rn , n,m > 2 and G u = Z 2 x Z 2 (G u = Z n , 
72 > 2 or Z n x Z m , 72 , m > 2) for Majorana (Dirac) neutrinos, the matrices C/ e and U u are hxed 
(up to diagonal phase matrices on the right, which are either unphysical for Dirac neutrinos, 
or contribute to the Majorana phases otherwise, and permutations of columns) by the residual 
symmetries of the charged lepton and neutrino mass matrices. In the case when the discrete 
symmetry Gf is fully broken in one of the two sectors, the corresponding mixing matrix U e 
or U u is unconstrained and contains in general three angles and six phases. 

Our article is organised as follows. In Section [2] we describe the parametrisations of 
the PMNS matrix depending on the residual symmetries G e and G v considered above. In 
Sections 10 and [5] we consider the breaking patterns 1), 2), 3) and derive sum rules for 
cos S. At the end of each of these sections we present numerical predictions for cos <5 in the 
cases of the flavour symmetry groups Gf = A 4 , T' , £4 and A 5 . In Section [b] we provide a 
summary of the sum rules derived in Sections i-m Further, in Sections [7] and [8] we derive 
the sum rules for the cases 4) and 5), respectively. In these cases the value of cos 6 cannot 
be fixed without additional assumptions on the unconstrained matrix U e or U v . The cases 
studied in 14 belong to the ones considered in Section [7J where the particular forms of the 
matrix U e , leading to sum rules of interest, have been considered. In Section[9]we present the 
summary of the numerical results. Section 10 contains the conclusions. Appendices [~\| [H| [C| 
[D] and [E] contain technical details related to the study. 


2 Preliminary Considerations 


As was already mentioned in the Introduction, the residual symmetries of the charged lepton 
and neutrino mass matrices constrain the forms of the matrices U e and U v and, thus, the form 
of the PMNS matrix U. To be more specific, if the charged lepton mass term is written in the 
left-right convention, the matrix U e diagonalises the hermitian matrix M e Me, UlM e MlU e = 
diag(m.g, mjp m^), M e being the charged lepton mass matrix. If G e is the residual symmetry 

group of M e Ml we have: 

p{g e )^M e Mlp{g e ) = M e M\ , (13) 


where g e is an element of G e , p is a unitary representation of Gf and p(g e ) gives the action 
of G e on the LH components of the charged lepton fields having as mass matrix M e . As 
can be seen from eq. (13), the matrices p(g e ) and M e Mj commute, implying that they are 
diagonalised by the same matrix U e . 

Similarly, if G u is the residual symmetry of the neutrino Majorana mass matrix M v one 
has: 

p(g u ) T M u p(g u ) = M v , (14) 


where g v is an element of G u , p is a unitary representation of Gf under which the LH flavour 
neutrino fields l = e, p, r, transform, and p(g u ) determines the action of G u on uuXx). 
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rt i\/f jsj commute, 

and therefore they can be diagonalised simultaneously by the same matrix U v . In the case of 


It is not difficult to show that also in this case the matrices p(g u ) and 


Dirac neutrinos eq. (14) is modified as follows: 


p{g v )'MlM„p(g v ) = M$M V 


(15) 


The types of residual symmetries allowed in this case and discussed below are the same as 
those of the charged lepton mass term. 

In many cases studied in the literature (e.g., in the cases of Gf = S 4 , A 4 , T', A§) p(gf), 
gf being an element of Gf, is assumed to be a 3-dimensional representation of Gf because 
one aims at unification of the three flavours (e.g., three lepton families) at high energy scales, 
where the flavour symmetry group Gf is unbroken. 

At low energies the flavour symmetry group Gf has necessarily to be broken to residual 
symmetries G e and G u , which act on the LH charged leptons and LH neutrinos as follows: 

II -t p(g e )lL , v lL ->• p{g v )v lL , 


where g e and g u are the elements of the residual symmetry groups G e and G v , respectively, 
and l L = (, e L ,fL L ,T L ) T , v lL = (v eL , v t l) T ■ 

The largest possible exact symmetry of the Majorana mass matrix M u having three non¬ 
zero and non-degenerate eigenvalues, is a Z 2 x Z 2 X Z 2 symmetry. The largest possible exact 
symmetry of the Dirac mass matrix M e is 17(1) x 17(1) X 17(1). Restricting ourselves to the 
case in which Gf is a subgroup of 517(3) instead of 17(3), the indicated largest possible exact 
symmetries reduce respectively to Z 2 x Z 2 and 17(1) x 17(1) because of the special determinant 
condition imposed from 51/(3). The residual symmetries G e and G u , being subgroups of Gf 
(unless there are accidental symmetries), should also be contained in 1/(1) x 1/(1) and Z 2 x Z 2 
(1/(1) x 1/(1)) for massive Majorana (Dirac) neutrinos, respectively. 

If G e = Z n , n > 2 or Z n x Z m , n,m > 2, the matrix U e is fixed by the matrix p(g e ) 
(up to multiplication by diagonal phase matrices on the right and permutations of columns), 
U e = U°. In the case of a smaller symmetry, i.e., G e = Z 2 , U e is defined up to a 1/(2) 
transformation in the degenerate subspace, because in this case p(g e ) has two degenerate 
eigenvalues. Therefore, 

u e = u 0 e u ij {e e ij ,5 e ij )^ k ^ u 

where Uij is a complex rotation in the i-j plane and \Ih are diagonal phase matrices, 

'hi = diag (e^ 1 ,1, l) , T 2 = diag (l, e^ 2 , l) , T 3 = diag (l, 1, . (16) 

The angle 6^- and the phases 5^, ipi, ip 2 and ^3 are free parameters. As an example of the 
explicit form of 17,we give the expression of the matrix Ui2(0f 2 , ^u)'- 

/ cos0“ 2 sin 0 “ 2 e -l<5 i 2 0\ 

^ 12 (^ 2 ^ 12 )= "Sin 0f 2 e w » cosdf 2 0 , (17) 

\ 0 0 1 / 

where a = e,v,o. The indices e, v indicate the free parameters, while “o” indicates the angles 
and phases which are fixed. The complex rotation matrices 1 / 23 (^ 23 ^ 23 ) an d ^ 13(^135 ^ 13 ) 

8 The right-left convention for the neutrino mass term is assumed. 
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are defined in an analogous way. The real rotation matrices Rij{6 ?•) can be obtained from 
Uij(6"j. 5fj) setting 5?- to zero, i.e., Rij(9 “•) = Uij(0fj, 0). In the absence of a residual symmetry 
no constraints are present for the mixing matrix U e , which can be in general expressed in terms 
of three rotation angles and six phases. 

Similar considerations apply to the neutrino sector. If G v = Z n , n > 2 or Z n x Z rn , 
n,m>2 for Dirac neutrinos, or G u = Z 2 x Z 2 for Majorana neutrinos, the matrix U v is fixed 
up to permutations of columns and right multiplication by diagonal phase matrices by the 
residual symmetry, i.e., U u = U°. If the symmetry is smaller, G v = Z 2 , then 

Obviously, in the absence of a residual symmetry, U u is unconstrained. In all the cases con¬ 
sidered above where G e and G u are non-trivial, the matrices p(g e ) and p{g u ) are diagonalised 
by U° and U°: 

(U°Jp(g e )U° e = p{g e ) d ^ and p{g„)U° v = p{g u ) d ^ . 

In what follows we define U° as the matrix fixed by the residual symmetries, which, in 
general, gets contributions from both the charged lepton and neutrino sectors, U° = (t/°)l[/°. 
Since U° is a unitary 3x3 matrix, we will parametrise it in terms of three angles and six 
phases. These, however, as we are going to explain, reduce effectively to three angles and one 
phase, since the other five phases contribute to the Majorana phases of the PMNS mixing 
matrix, unphysical charged lepton phases and/or to a redefinition of the free parameters 
contained in U e and U v . Furthermore, we will use the notation Of- , 6 F-, Sf-, <5'/ for the free 
angles and phases contained in U, while the parameters marked with a circle contained in 
U °, e.g., 0°j, 5°^, are hxed by the residual symmetries. 

In the case when G e = Z 2 and G u = Z n , n > 2 or Z n X Z m , n, m > 2 for massive Dirac 
neutrinos, or G v = Zo x Z 2 for Majorana neutrinos, we have: 

u = u v (et r 5 t^°u°(e^ o° l 3 , 02.3, {<%})Qo 
= s tj ~ ^)u°m 2 , 01 3 , e° 23 , {%})Q 0 , (18) 

where (ij) = ( 12 ), (13), (23) and {/>)/} = {^i 2) ^ 13 , ^ 23 }- The unitary matrix U° contains three 
angles and three phases, since the additional three phases can be absorbed by redefining the 
charged lepton fields and the free parameter <5®. (see below). Here is a diagonal matrix 
containing a hxed phase in the j-th position. Namely, 

*1 = diag (e^, 1, l) , ^2 = diag (l, e^, l) , T 3 = diag (l, 1, e^) . (19) 

The matrix Qo, defined in eq. ([ 3 ]), is a diagonal matrix containing two free parameters con¬ 
tributing to the Majorana phases. Since the presence of the phase 1 /)? amounts to a redefinition 
of the free parameter Sfj , we denote (Sfj — ip ?) as Sfj . This allows us to employ the following 
parametrisation for U: 


U = Uijietj, 8^)11°m 2 , ^ 13 ) 023, Wo, (20) 

where the unphysical phase matrix T? on the left has been removed by charged lepton re¬ 
phasing and the set of three phases {/))/} reduces to only one phase, 5^, since the other two 
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contribute to redefinitions of Qo, Sfj and to unphysical phases. The possible forms of the 
matrix U°, which we are going to employ, are given in Appendix |Bj 

For the breaking patterns G e = Z n , n > 2 or Z n x Z m , n,m > 2 and G u = Z 2 , valid for 
both Majorana and Dirac neutrinos, we have: 


= u°(9° 12 ,9° 13 ,9° 23 ,s° kl )u ij (9" j ,6” j - 


&ii)Qo 

^ 0 + ^°)^°Qo, 


( 21 ) 


where (ij ) = (12), (13), (23), and the two free phases, which contribute to the Majorana phases 
of the PMNS matrix if the massive neutrinos are Majorana particles, have been included in 
the diagonal phase matrix Qo- Notice that if neutrinos are assumed to be Dirac instead 
of Majorana, then the matrix Qo can be removed through re-phasing of the Dirac neutrino 
fields. Without loss of generality we can redefine the combination 5 T — ij)° + -0? as <$- and 
the combination T°T°Qo as Qo, so that the following paranretrisation of U is obtained: 

U = U°(9° 12 ,91 3 ,0 2 ° 3 , ShWijWj^Qo . (22) 


In the case of G e = Z 2 and G v = Z 2 for both Dirac and Majorana neutrinos, we can write 


U = U ij (9^dt j )^ j U 0 (9° 12 , 91 3 ,0 2 ° 3 ,8 ° kl )*° r *°U rs W s , 5" S )Q 0 

= vpijie&sfj - r j )u°(9° v 2 , o° 13 ,e° 23 , s° kl )u rs (e , v a ,c - w° + r a W r KQo 


(23) 


with (ij) = (12), (13), (23), (rs) = (12), (13), (23). The phase matrices T° are defined as in 
eq. ( |19[ ). Similarly to the previous cases, we can redefine the parameters in such a way that 
U can be cast in the following form: 


U = Uiiietj, Sfj)U 0 (6l 2 ,91 3 , 91 3 , 5° kl )U rs (9^ OQo , (24) 

where Qo can be phased away if neutrinos are assumed to be Dirac particles j^] 

If G e is fully broken and G u = Z n , n > 2 or Z n x Z m , n, m > 2 for Dirac neutrinos or 
G v = Z 2 x Z 2 for Majorana neutrinos, the form of U reads 


U = U(9l 2l 9\ 3 , 91 3 , 5 e rs )W 3 U°(9° 12l 9° 13 , 0 2 ° 3 , {6° kl })Q 0 , (25) 

Notice that in general we 


(16) 


where the phase matrices T 2 and T 3 are defined as in eq. 
can effectively parametrise U° in terms of three angles and one phase since of the set of three 
phases two contribute to a redefinition of the matrices Qo, ^2 and ^ 3 . Furthermore, 

under the additional assumptions on the form of U(9f 2 , 9 f 3 , 9 23 , S^ s ) and also taking {6 kl } = 0, 

In the numerical 


analyses performed in 

11 

If 

14 , the angles 9°- have been set, in particular, to the values 

corresponding to the T1 

BM, B] 

vl (LC), GRA, GRB and HG symmetry forms. 


Finally for the breaking patterns G e = Z n , n > 2 or Z n x Z m , n,m > 2 and G u fully 
broken when considering both Dirac and Majorana neutrino possibilities, the form of U can 


be derived from eq. (25) by interchanging the fixed and the free parameters. Namely, 


The cases found in 


U = U°(9° 12 ,9° 13 ,9° 23 , S° kl )W 3 U(9 v 12 ,9 
eqs. pOb, p|, d24b, p5| and p| 


13) ^23) Ks)Q 0 • 
are summarised in 


(26) 

Table [l| The 


reduction of the number of free parameters indicated with arrows corresponds to a redefinition 
of the charged lepton fields. 


9 We will not repeat this statement further, but it should be always understood that if the massive neutrinos 
are Dirac fermions, then two phases in the matrix Qo are unphysical and can be removed from U by a re-phasing 
of the Dirac neutrino fields. 
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G e cG f 


fully broken 

Z 2 

Z n , n > 2 or 
Z n x Z m , n,m> 2 

fully broken 
fully broken 

Z‘2 

Z n , n > 2 or 
Z n x Z m , n,m> 2 

Z 2 


Z n , n > 2 or 
x Z m , n, m > 2 


G u cG f 

fully broken 
fully broken 

fully broken 

Z 2 

Z n . n > 2 or 
Z n x Z m , n, m > 2 

Z 2 

Z 2 


U e d.o.f. 

U v d.o.f. 

U d.o.f. 

9—^6 

9 —8 

12 ^ 4 (+2) 

4—^2 

9 —8 

10 ^ 4 (+2) 

0 

9 —>■ 8 

8-^4 (+2) 

9-^6 

4 

10 -> 4 (+2) 

9—^6 

2 

8-^4 (+2) 

4—^2 

4 

4 (+2) 

0 

4 

2 (+2) 


Z n , n > 2 or 
Z n x , n,m >2 

Z n . n > 2 or 
Z n x Z m , n, m > 2 


2 (+ 2 ) 
0 (+ 2 ) 


Table 1: Number of effective free parameters, degrees of freedom (d.o.f.), contained in U 
relevant for the PMNS angles and the Dirac phase (and Majorana phases) in the cases of the 
different breaking patterns of Gf to G e and G v . Arrows indicate the reduction of the number 
of parameters, which can be absorbed with a redefinition of the charged lepton fields. 


In the breaking patterns considered, it may be also possible to impose a generalised CP 
(GCP) symmetry. An example of how imposing a GCP affects the sum rules is shown in 
Appendix [Dj In the case in which a GCP symmetry is preserved in the neutrino sector we 
have for the neutrino Majorana mass matrix 


39 


X M V X { = Ml 


(27) 


Since the matrix X % is symmetric there exists a unitary matrix such that Xj = Qf and 
QjM v i\ is real. Therefore when GCP is preserved in the neutrino sector, the phases in the 
matrix U v can be fixed. An alternative possibility is that GCP is preserved in the charged 
lepton sector, which leads to the condition 


39 


(A'f)t M e M\X\ = {M e M\y 


(28) 


Since ( Xf) T = Xf, the phases in the matrix U e are fixed, because (Qf)tM e MjQ? is real. 
The fact that the matrices Xi, if GCP is preserved in the neutrino sector, or Xf if it is 
preserved in the charged lepton sector, are symmetric matrices can be proved applying the 


GCP transformation twice. In the first case, eq. (27) allows one to derive the general form of 


9 








Xi 40-42 : 


while in the latter case 


Xi = u v xf ia *uT. 


Xf = U e (X^Uj. 


(29) 


(30) 


Equations (29) and (30) imply that Xi and Xf are symmetric matrices 


10 


We note finally that the titles of the following sections refer to the residual symmetries of 
the charged lepton and neutrino mass matrices, while the titles of the subsections reflect the 


free complex rotations contained in the corresponding parametrisation of U, eqs. (20), (22), 


(24), (25) and (26). 


3 The Pattern G e = Zo and G v = n > 2 or Z n x Z m , n, m > 2 


In this section we derive sum rules for cos 5 for the cases given in eq. (20). Recall that the 


matrix U e is fixed up to a complex rotation in one plane by the residual G e = Z 2 symmetry, 
while U v is completely determined (up to multiplication by diagonal phase matrices on the 
right and permutations of columns) by the G v = Z 2 y Z 2 residual symmetry in the case 
of neutrino Majorana mass term, or by G u = Z n , n > 2 or Z n x Z m , n,m > 2, residual 
symmetries if the massive neutrinos are Dirac particles. At the end of this section we will 
present results of a study of the possibility of reproducing the observed values of the lepton 
mixing parameters sin 2 6 12 , sin 2 $13 and sin 2 $23 and of obtaining physically viable predictions 
for cos 5 in the cases when the residual symmetries G e = Z 2 and G u = Z n , n > 2 or Z n x Z m , 
n,m> 2, originate from the breaking of the lepton flavour symmetries A 4 (T 7 ), S 4 and A 5 . 

3.1 The Case with ^12(^12^12) Complex Rotation (Case Al) 


Employing the parametrisation of the PMNS matrix U given in eq. (20) with (ij) = (12) and 
the parametrisation of U° given as 


($121 ^ 13 > ^ 23 ! ^12) — ^12(^12)^12)^23(^23)^13(^13)’ 
we get for U (see Appendix [B| for details): 


(32) 


u = u 1 2(9t2Zt2)Ui2(ei2,s°i2)R2 3 (0 o 23 )R 13 (e o 13 )Q 0 . (33) 

The results derived in Appendix [B] and given in eq. ( | 212 ) allow us to cast eq. (33) in the form: 

u = Ri2(di2)Pi(6 1 2)R23(0°2 3 )Ri 3 (6° 13 )Q 0 , Pi(5i 2 ) = diag(e^ 12 ,1,1), (34) 


with < 5 i 2 = a — f3, where sin$ 12 , a and /3 are defined as in eqs. (213) and (214) after setting 
i = 1, j = 2, 0f 2 = 0 i 2 ) $12 = $® 2 ) 0 i 2 = 012 an d $12 = $ 12 ' Using eq. (34) and the standard 


10 This fact can be also derived from the requirement that the GCP transformations contain the physical CP 
transformation, i.e., the GCP transformations applied twice to a field should give the field itself [40 43[44| : 

4>{x) —1 X T (j>*(xp) X r X*(f>(x) = 4>{x), (31) 

where x = (xo,x), x v = (xo, —x). The notation we have used for X r emphasises the representation r for the 
GCP transformations. 


10 






























parametrisation of the PMNS matrix U, we find: 


sin 2 0 i 3 = | C/ e 31 2 = cos 2 0 i 2 sin 2 0 f 3 + cos 2 0 f 3 sin 2 0 i 2 sin 2 0 ^ 


+ - sin 2012 sin 20 f 3 sin 9^ cos hi 2 , 


sin~ 0 23 = 


sin^ 0i2 = 


I U, 


M3| 


1 


1 - \u e3 \ 2 
i - \u e3 \ 2 


COS 2 013 

-,2 


[ sin 2 0 i 3 — sin 2 0 33 + cos 2 0 f 3 sin 2 0 ^] , 


cos 2 023 sin 2 012 

COS 2 013 


(35) 

(36) 

(37) 


From eqs. (35) and (36) we get the following correlation between the values of sin 2 0i 3 and 
sin 2 0 23 : 

sin 2 0 i 3 + cos 2 0 i 3 sin 2 02 3 = sin 2 0f 3 + cos 2 0f 3 sin 2 9^ . (38) 


Notice that eq. (37) implies that 


sin 2 0 i 2 = 


cos 2 0 i 3 sin 2 0 i 2 


cos 


2 0 , 


(39) 


23 


In order to obtain a sum rule for cos h, we compare the expressions for the absolute 
value of the element U T 2 of the PMNS matrix in the standard parametrisation and in the 
parametrisation defined in eq. (34), 


| U T 21 = | cos 0i2 sin 023 + sin 0i 3 cos 02 3 sin 9\ 2 e l ° | = | sin 0^ 3 1 . 

From the above equation we get for cosh: 

cos 2 0i 3 (sin 2 0^3 — cos 2 0i2) + cos 2 0f 3 cos 2 02 3 (cos 2 0i2 — sin 2 0i2 sin 2 0i 3 ) 


cosh = 


sin 2012 sin 0i 3 1 cos 0f 3 cos 0o 3 1 (cos 2 0i 3 — cos 2 0f 3 cos 2 02 3 ) 2 


(40) 


(41) 


For the considered specific residual symmetries G e and G u , the predicted value of cosh in the 
case A1 discussed in this subsection depends on the chosen discrete flavour symmetry Gf via 
the values of the angles 0 f 3 and 02 3 . 

The method of derivation of the sum rule for cosh of interest employed in the present 
subsection and consisting, in particular, of choosing adequate parametrisations of the PMNS 
matrix U (in terms of the complex rotations of U e and of U u ) and of the matrix U° (determined 
by the symmetries G e , G v and Gf), which allows to express the PMNS matrix U in terms of 
minimal numbers of angle and phase parameters, will be used also in all subsequent sections. 
The technical details related to the method are given in Appendices [B] and [C] 

We note finally that in the case of hf 2 = 0, the symmetry forms TBM, BM, GRA, GRB 
and HG can be obtained from U° = -Ri 2(0 i 2)7?23(023)77 l 3(0 i 3 ) for specific values of the angles 
given in Table [2| In this case, the angles 9°■ are related to the angles 9 V VJ defined in Section 2.1 


of ref. 14 as follows: 


• 2 /io 2 r\v • 2 

sm #23 = cos # 12 sin # 2 3 ? 


sin 


7 13 


sin 2 023 sin 2 0^ 2 


1 — sin 2 9, 


sin 2 0 f 2 = 


sin 


T2 


23 


1 — sim 


(42) 


23 
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Mixing 

0°V2 

0 23 

Oh 

TBM 

7t/4 

— sin- 1 (l/v / 3) 

n/6 

BM 

sin- 1 ^2/3 

—7t/6 

sin” 1 (l/y/3) 

GRA 

sin- 1 yj( 7- n/5)/11 

-sin’ 1 ^(5 + y/5)/20 

sin" 1 s]{7 - x/5)/22 

GRB sin' 1 ^2(15 - 2y/5)/41 

— sin- 1 ^/(3 + \/5) /16 

sin” 1 y/(15 - 2>/5)/41 

HG 

sin 1 ^2/5 

— sin -1 ^3/8 

sin -1 i/l/5 


Table 2: The symmetry forms TBM, BM (LC), GRA, GRB and HG obtained in terms of the 
three rotations -Ri 2 ( 0 i 2 )R 23 ( 02 3 )Ri 3 ( 0 ° 3 )- 


3.2 The Case with U\ 3 {6l 3 ,5l 3 ) Complex Rotation (Case A2) 


Using the parametrisation of the PMNS matrix U given in eq. 
following parametrisation of U°, 


u°{e \ 2 , e° 13 , o° 23l s° 13 ) = u 13 (9° 13 , 5° 13 )R 23 (e° 23 )R 12 (e° 12 ), 

we get for U (for details see Appendix [B]) : 

U = C 13 (^3,^3)^13(^3,^3)^23(0 2 O 3 ) J Rl2(0^)Qo - 


with ( ij ) = (13) and the 
(43) 


(44) 


The results derived in Appendix [B| and presented in eq. (212) allow us to recast eq. (44) in 
the following form: 


u = i?i 3 (0i 3 )Pi(^i 3 ) J R23(023)^i2(0i O 2)Qo , Pi(6ia) = diag(e^ 13 ,1,1). 


(45) 


Here <h 3 = 0.-/3, where sin0i 3 , a and /3 are defined as in eqs. (213) and (214) after setting 
* = 1, j = 3, 0“ 3 = 0f 3 , <5“ 3 = <5f 3 , 9 33 = 0° 3 and S 33 = d^ 3 . Using eq. (45) and the standard 
parametrisation of the PMNS matrix U, we find: 


sin 2 0i 3 = \U e3 \= sin 2 0 i3 cos 2 02 3 , 


sin 2 0 23 = 


I ^31 2 


sin 2 0° 3 


sin 2 0 12 = 


l-|C e 3| 2 1-sin 2 013 ’ 

|C e2 | 2 1 


l-|C e3 | 2 1-sin 2 0 i 3 


[ cos 2 0i3 sin 2 0^2 + cos 2 0^ 2 sin 2 0 33 sin 2 9 23 


1 


— - sin 2013 sin 20\ 2 sin 0^ 3 cos di 3 ] . 


(46) 

(47) 

(48) 


Thus, in this scheme, as it follows from eq. (47), the value of sin" 0 23 is predicted once the 
symmetry group Gf is fixed. This prediction, when confronted with the measured value of 
sin 2 023 , constitutes an important test of the scheme considered for any given discrete (lepton 
flavour) symmetry group Gf, which contains the residual symmetry groups G e = Z 2 and 
G u = Z n , n > 2 and/or Z n x Z m , n,m> 2 as subgroups. 
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As can be easily demonstrated, the case under discussion coincides with the one analysed 
in Section 2.2 of ref. |14j. The parameters 0 23 and 9 \ 2 in [14j can be identified with $33 and 0f 2 , 
respectively. Therefore the sum rule we obtain coincides with that given in eq. (32) in |14|: 


cos 5 = 


cos 2 0 i 3 (cos 2 9 


° 2 cos 2 023 — cos2 012 ) + sin 2 023 (cos 2 0 i 2 — sin 2 0 i 2 sin 2 013 ) 
sin 2012 sin 0 i 31 sin 0 23 |(cos 2 0 i 3 — sin 2 0 ^ 3 ) 2 


(49) 


The dependence of cos <5 on Gf in this case is via the values of the angles 0 ° 2 and 9 23 


3.3 The Case with U 23 [9 23 , 5 23 ) Complex Rotation (Case A3) 

In the case with ( ij ) = (23), as can be shown, cos<5 does not satisfy a sum rule, i.e., it 
cannot be expressed in terms of the three neutrino mixing angles 0 i 2 , 0 i 3 and 023 and the 
other fixed angle parameters of the scheme. Indeed, employing the parametrisation of U° as 
t/°( 0 ? 2 , 0i 3 , 0 2 3 , ^ 3 ) = ^ 23 ( 023 ,^ 23 )^ 13 ( 013 )^ 12 ( 012 )) we can write the PMNS matrix in the 
following form: 


U = U 23 (9 e 23 ,8 e 23 )U 23 {9° 23 , 5° 23 )R 13 (9° 13 )R 12 (9° 12 )Q 0 . (50) 

Using the results derived in Appendix |B| and shown in eq. (212), we can recast eq. (50) as 
U = R 2 3 (. 9 23 )P 2 ( 5 23 )R 13 ( 9 ° 13 )R 12 ( 9 ° 12 )Q 0 , P 2 (S 23 ) = diag(l, e^ 23 ,1), (51) 


with 5 23 = a — f3 , where sin 023, a and /3 are defined as in eqs. (213) and (214) after setting 
i = 2 , 3 = 3, 0f 3 = 0| 3 , #? 3 = 5| 3 , 9 b 2 3 = 0 23 and S% 3 = S 23 . Comparing eq. ([m} and the 
standard parametrisation of the PMNS matrix, we find that sin 2 013 = sin 2 0 J 3 , sin 2 0 23 = 
sin 2 0 23 , sin 2 0 i 2 = sin 2 9 32 and cos 8 = ± cos 5 23 . 

It follows from the preceding equations, in particular, that since, for any given Gf com¬ 
patible with the considered residual symmetries, 0^3 and 0^2 have fixed values, the values of 
both sin 2 0 i 3 and sin 2 0 i 2 are predicted. The predictions depend on the chosen symmetry Gf. 
Due to these predictions the scheme under discussion can be tested for any given discrete 
symmetry candidate Gf, compatible, in particular, with the considered residual symmetries. 

We have also seen that 5 is related only to an unconstrained phase parameter of the 
scheme. In the case of a flavour symmetry Gf which, in particular, allows to reproduce 
correctly the observed values of sin 2 0 i 2 and sin 2 0 i 3 , it might be possible to obtain physically 
viable prediction for cos <5 by employing a GCP invariance constraint. An example of the effect 
that GCP invariance has on restricting CPV phases is given in Appendix [D} Investigating the 
implications of the GCP invariance constraint in the charged lepton or the neutrino sector in 
the cases considered by us is, however, beyond the scope of the present study. 

3.4 Results in the Cases of Gf = A 4 ( T'), S 4 and A 5 


The cases detailed in Sections 3.1 - 3.3 can all be obtained from the groups A 4 (T '), S 4 and 


A 3 , when breaking them to G e = Z 2 and G v = Z n (n > 3) in the case of Dirac neutrinos, or 
G u = Z 2 x Z 2 in the case of both Dirac and Majorana neutrinos^] We now give an explicit 
example of how these cases can occur in A 4 . 

In the case of the group A 4 (see, e.g., [45]), the structure of the breaking patterns discussed, 
e.g., in subsection 3.1 can be realised when i) the S generator of A 4 is preserved in the 


1 We only consider Z2 x Z2 when it is an actual subgroup of Gf. 
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neutrino sector, and when, due to an accidental symmetry, the mixing matrix is fixed to be 
tri-bimaximal, U° = C/tbMj up to permutations of the columns, and ii) a Z 2 2st or Z 2 ^ Tl 
is preserved in the charged lepton sector. The group element generating the Z -2 symmetry is 
diagonalised by the matrix U°. Therefore the angles 9l 2 , $13 and 9 23 are obtained from the 
product U° = (U°)^U°. The same structure (the structure discussed in subsection 3.2) can 
be obtained in a similar manner from the flavour groups S4 and A3 {A4, S4 and A5). 

We have investigated the possibility of reproducing the observed values of the lepton mix¬ 
ing parameters sin 2 $ 12 , sin 2 $13 and sin 2 023 as well as obtaining physically viable predictions 
for cos 5 in the cases of residual symmetries G e = Z2 and G v = Z n , n > 2 or Z n x Z m , 
n,m > 2 12 (Dirac neutrinos), or G u = Z2 z Z2 (Majorana neutrinos), discussed in sub¬ 


sections and |3.3| denoted further as Al, A2 and A3, assuming that these residual 

symmetries originate from the breaking of the flavour symmetries A4 ( T '), S4 and A3. The 
analysis was performed using the current best fit values of the three lepton mixing parameters 
sin 2 0 i 2 , sin 2 0 i 3 and sin 2 023. The results we have obtained for the symmetries A4 ( T'), S4 
and A 3 are summarised below. 

We have found that in the cases under discussion, i.e., in the cases Al, A2 and A3, and 
flavour symmetries Gf = A4 ( T'), S4 and A3, with the exceptions to be discussed below, 
it is impossible either to reproduce at least one of the measured values of sin 2 0 i 2 , sin 2 0 i 3 
and sin 2 023 even taking into account its respective 3a uncertainty, or to get physically viable 
values of cos 5 satisfying |cos<5| < 1. In the cases Al and A2 and the flavour groups A 4 
and S4, for instance, the values of cos <5 are unphysical. Using the group Gf = A3 leads 
either to unphysical values of cos <5, or to values of sin 2 023 which lie outside the corresponding 


current 3tr allowed interval. In the case A3 (discussed in subsection 3.3), the symmetry A4, 
for example, leads to (sin 2 012 , sin 2 0i 3 ) = (0,0) or (1,0). 

As mentioned earlier, there are three exceptions in which we can still get phenomenologi¬ 
cally viable results. In the Al case (A2 case) and S 4 flavour symmetry, one obtains bimaximal 
mixing corrected by a complex rotation in the 1-2 plane p*| (1-3 plane). The PMNS angle 023 
is predicted to have a value corresponding to sin 2 623 = 0.488 (sin 2 623 = 0.512). For the best 
fit values of sin 2 0i2 and sin 2 0i3 we find that cos 5 = — 1.29 (cos 6 = + 1.29). However, using 
the value of sin 2 0i2 = 0.348, which lies in the 3d allowed interval, one gets the same value 
of sin 2 023 and cos <5 = —0.993 (cos 5 = 0.993), while in the part of the 3cr allowed interval of 
sin 2 0i2, 0.348 < sin 2 0i2 < 0.359, we have —0.993 < cos 5 < —0.915 (0.993 > cos 6 > 0.915). 

Also in the Al case (A2 case) but with an A3 flavour symmetry and residual symmetry 
G u = Z 3 , which is only possible if the massive neutrinos are Dirac particles, we get the 
predictions sin 2 023 = 0.553 (sin 2 023 = 0.447) and cos 5 = 0.716 (cos <5 = —0.716). In the Al 
case (A2 case) with an A 3 flavour symmetry and residual symmetry G u = Z 3 , which can be 
realised for neutrino Dirac mass term only, for the best fit values of sin 2 0i2 and sin 2 0i3 we 
get the predictions sin 2 623 = 0.630 (sin 2 623 = 0.370), which is slightly outside the current 3er 
range) and cos <5 = — 1.12 (cos <5 = 1.12). However, using the value of sin 2 0i2 = 0.321, which 
lies in the la allowed interval of sin 2 0i2, one gets the same value of sin 2 023 and cos 5 = —0.992 


12 Note that there are no subgroups of the type Z n x Z m bigger than Z 2 x Z 2 in the cases of A 4 , S 4 and A 5 . 
13 For the case Al it can been shown that 


diag(—1,1, 1)U(0° 12 ,6° 12 )R(0° 23 )R(0° 13 ) diag(l, - 1 , 1 ) = Ubm , 


(52) 


if $23 = sin -1 (1 /2), #13 = sin - 1 (-y/l/3), 0 ° 2 = tan -1 (-y/3/2 + 1/2) and 842 = 0. Therefore, one has BM 

mixing corrected from the left by a 1/(2) transformation in the degenerate subspace in the 1-2 plane. Note 
that our results are in agreement with those obtained in 46 . 


14 









(cos 5 = 0.992). In the part of the 3 a allowed interval of sin 2 0 i 2 , 0.321 < sin 2 6*12 < 0.359, 
one has —0.992 < cos 6 < —0.633 (0.992 > cosd > 0.633). 


4 The Pattern G e = Z n , n > 2 or Z n x Z m , n, m > 2 and G u = Z 2 


In this section we derive sum rules for cos 6 in the case given in eq. (|22j). We recall that for 
G e = Z n , n > 2 or Z n x Z m , n,m> 2 and = Z 2 of interest, the matrix U e is unambiguously 
determined (up to multiplication by diagonal phase matrices on the right and permutations 
of columns), while the matrix U u is determined up to a complex rotation in one plane. 

4.1 The Case with ^13(^13,^13) Complex Rotation (Case Bl) 


Combining the parametrisation of the PMNS matrix U given in eq. (22) with (ij) = (13) and 
the parametrisation of U° as 

U°m 2 , 91 3 , 0 2 ° 3 , 61 3 ) = R 23 (9° 23 )R 12 (9° 12 )U 13 (9° 13 , S ° 13 ), (53) 

we get for U (the details are given again in Appendix [B| : 

U = • (54) 

The results derived in Appendix [B] and reported in eq. ( |212 ) allow us to recast eq. (54) in the 
form: 

U = R23(0° 23 )Ri2(0 o 12 )P 3 (5 13 )R 13 (d 13 )Q 0 , P 3 (6 13 ) = diag(l, 1, e^ 13 ). (55) 

Here <$13 = — a —/3 and we have redefined ^ 13 ( 0 , /3)Qo as Qo, where p 3 (a , (3) = diag(e*“, 1 , e^ 3 ) 
and the expressions for sin 2 ^ 13 , a and (3 can be obtained from eqs. (213) and (2141), by setting 
i = 1, j = 3, &i 3 = 91 3 , df 3 = < 5 ^ 3 , 3 = 9\ 3 and 5± 3 = 6^. Using eq. (55) and the standard 

parametrisation of the PMNS matrix U, we find: 

sin 2 0 i 3 = \U e3 \ 2 = cos 2 0i 2 sin 2 0 i3 , (56) 

' T Gs\ 2 _ 1 


sin 023 = 


l-\U e3 \ 2 COS 2 013 


[ cos 2 023 sin 2 0 i 3 sin 2 0 i 2 + cos 2 0 i 3 sin 2 0 


2 QO 
23 


1 


— - sin 2013 sin 20^ 3 sin 0 ^ cos 


< 5 1,3] 


sin 2 0 i 2 = 


We2\ 2 


sin 2 0 i 2 


l-\U e3 \ 2 COS 2 013 


(57) 

(58) 


It follows from eq. (58) that in the case under discussion the values of sin 2 0 i 2 and sin 2 0 i 3 are 
correlated. 

A sum rule for cos 6 can be derived by comparing the expressions for the absolute value of 
the element U t2 of the PMNS matrix in the standard parametrisation and in the one obtained 
using eq. (55): 

| C/t -2 I = | cos 0 i 2 sin 023 + sin 0 i 3 cos 0 23 sin 0 i 2 e ?5 | = | cos 0 J 2 sin 0^1 . 

From this equation we get 

cos 2 0 i 3 (cos 2 0 i 2 cos 2 023 — cos 2 0 23 ) + sin 2 0 i 2 (cos 2 6 23 — sin 2 0 i 3 sin 2 0 23 ) 
sin 2023 sin 0 i 31 sin 0 i 2 1 (cos 2 0 i 3 — sin 2 0 ° 2 ) 2 
The dependence of the predictions for cos <5 on Gf is in this case via the values of 0 ° 2 and 0 23 . 


(59) 


(60) 
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4.2 The Case with U 23 (9 23 , 5 23 ) Complex Rotation (Case B2) 

Utilising the parametrisation of the PMNS matrix U given in eq. \22\ with (ij) = (23) and 
the following parametrisation of U°, 


U°(9l 2 , ^ 3 , 01 3, <5 2 ° 3 ) = Ris(ei 3 )R 12 (e° 12 )U 23 (9° 23 ,5 ° 23 ), (61) 

we obtain for U (Appendix [B| contains the relevant details): 

U = Ri3(9° 13 )R 12 (9° 12 )U 23 (9° 23 ,6° 23 )U 23 (9 % 3 , S^)Q 0 . (62) 

The results given in eq. (212) in Appendix [B| make it possible to bring eq. (62) to the form: 

U = Ri3(9l 3 )R 12 (9° 12 )P 3 (S 23 )R 23 (9 23 )Q 0 , P 3 (5 23 ) = diag(l, 1, e ^). (63) 

Here < 5 2 3 = —a—/3 and we have redefined P 2 3 (a, /3)Qo as Qo, where U 23 (a, /3) = diag(l, e* a , e 1 ^). 
Using eq. (63) and the standard parametrisation of the PMNS matrix U, we find: 

sin 2 0 i 3 = \U e3 \ 2 = cos 2 0 i 3 sin 2 0 ° 2 sin 2 0 23 + sin 2 0i 3 cos 2 0 2 


723 


+ - sin 20 2 3 sin 20^ sin 9 32 cos < 5 2 3 , 


sin 2 023 = 
sin 2 0 i 2 = 




cos 2 0 i 2 sin 2 0 


23 


l-\U e3 \ 2 COS 2 013 

I U e2 | 2 _ COS 2 013 — COS 2 012 COS 2 013 

1 - |C e3 | 2 COS 2 013 


(64) 

(65) 

( 66 ) 


Equation (66) implies that, as in the case investigated in the preceding subsection, the values 
of sin 2 0 i 2 and sin 2 d 33 are correlated. 

The sum rule for cos <5 of interest can be obtained by comparing the expressions for the 
absolute value of the element U T \ of the PMNS matrix in the standard parametrisation and 
in the one obtained using eq. (63): 

| U T 1 1 = | sin 0i2 sin 9 23 — sin 0\ 3 cos 0i 2 cos 0 23 e lS \ = | cos 0J 2 sin 0i 3 1 . 

From the above equation we get for cos 5: 

cos 2 0 i 3 (sin 2 0 32 — cos 2 0 23 ) + cos 2 0 32 cos 2 0 33 (cos 2 023 — sin 2 9\ 3 sin 2 0 23 ) 
sin 2023 sin 0i 3 1 cos 0\ 2 cos 9 \ 3 1 (cos 2 0i 3 — cos 2 0 32 cos 2 0 33 ) 2 
The dependence of cos 5 on Gf is realised in this case through the values of 0 32 and 0i 3 


(67) 


( 68 ) 


4.3 The Case with Ui 2 {6 32 ,5 32 ) Complex Rotation (Case B3) 

In this case, as we show below, cos 5 does not satisfy a sum rule, and thus is, in gen¬ 
eral, a free parameter. Indeed, using the parametrisation of U° as t/°(0i 2 , 0J 3 , 0 23 , 5 32 ) = 
R 2 3(9 23 )Ri 3 (91 3 )Ui 2 (9 32 ,5 i 2 ) we get the following expression for U: 


U = R 2 3(9h)Ri3(9l 3 )U 12 (9l 2 ,6l 2 )Un(9 u 12 ,5r 2 )Q 0 . 


After recasting eq. (69) in the form 


U = R 2 3(9 0 23 )R 13 (9° 13 )P 2 (S 12 )R 12 (9 12 )Qo , P 2 (S 12 ) = diag(l, e^, 1), 


(69) 


(70) 
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where 812 = —a — /3, we find that sin 2 $13 = sin 2 0J 3 , sin 2 023 = sin 2 $231 sin 2 $12 = sin 2 0 i 2 and 

COS 5 = ± COS § 12 - 

It follows from the expressions for the neutrino mixing parameters thus derived that, given 
a discrete symmetry Gf which can lead to the considered breaking patterns, the values of 
sin 2 0 i 3 and sin 2 023 are predicted. This, in turn, allows to test the phenomenological viability 
of the scheme under discussion for any appropriately chosen discrete lepton flavour symmetry 


Gf 


In what concerns the phase 8 , it is expressed in terms of an unconstrained phase parameter 

is 


present in the scheme we are considering. The comment made at the end of subsection 3.3 


valid also in this case. Namely, given a non-Abelian discrete flavour symmetry Gf which allows 
one to reproduce correctly the observed values of sin 2 0 i 3 and sin 2 O 23 , it might be possible to 
obtain physically viable prediction for cos 8 by employing a GCP invariance constraint in the 
charged lepton or the neutrino sector. 


4.4 Results in the Cases of Gf = A 4 (T'), S 4 and A 5 


The schemes discussed in Sections 4.1 4.3 are realised when breaking Gf = A 4 ( T £4 and 

A3, to G e = Z n (n > 3) or Z 2 x Z 2 and G v = Zo, for both Dirac and Majorana neutrinos. As 
a reminder to the reader, we investigate the case of Z 2 x Z 2 when it is an actual subgroup 
of Gf. As an explicit example of how this breaking can occur, we will consider the case of 
Gf = A4 ( T '). The other cases when Gf = £4 or A3 can be obtained from the breaking of £4 
and A 3 to the relevant subgroups as given in 46 and 47 , respectively. 

In the case of the group A 4 (see, e.g., |45|), the structure of the breaking patterns discussed, 


e.g., in subsection 4.1 can be obtained by breaking A 4 i) in the charged lepton sector to any of 


the four Z 3 subgroups, namely, Zjf, Zf T , Zj s , Z ^ TS , and ii) to any of the three Z 2 subgroups, 
namely, Z:f, Zj~ ST , Z ^ st2 , in the neutrino sector. In this case the matrix U° = I/tbm gets 
corrected by a complex rotation matrix in the 1-3 plane coming from the neutrino sector. 

The results of the study performed by us of the phenomenological viability of the schemes 
with residual symmetries G e = Z n , n > 2 or Z n x Z m , n,m> 2 and G u = Z 2 , discussed in 


subsections 4.1 4.2 and 4.3, and denoted further as Bl, B2 and B3, when the residual sym¬ 
metries result from the breaking of the flavour symmetries A 4 (T'), £4 and A 3 , are described 
below. We present results only in the cases in which we obtain values of sin 2 0i2, sin 2 0i3 
and sin 2 623 compatible with their respective measured values (including the corresponding 
3cr uncertainties) and physically acceptable values of cos 8 . 

For Gf = A4, we find that only the case Bl with G e = Z3 is phenomenologically viable. 
In this case we have (sin 2 0 32 ,sin 2 0 23 ) = (1/3,1/2), which leads to the predictions sin 2 0 i 2 = 
0.341 and cos <5 = 0.570. We find precisely the same results in the case Bl if Gf = £4 and 
G e = Z3. Phenomenologically viable results are obtained for Gf = £4 and G e = Z3 in 
the case B2 as well. In this case (sin 2 0^ 2 , sin 2 0J 3 ) = (1/6,1/5), implying the predictions 
sin 2 0 i 2 = 0.317 and cos<5 = —0.269. If G e = Z4 or Z2 x Z2 results from Gf = £4, we get 
in the case Bl (sin 2 0 32 , sin 2 0 23 ) = (1/4,1/3) and correspondingly sin 2 0 i 2 = 0.256 (which 
lies slightly outside the current 3 a allowed range of sin 2 012 ) and the unphysical value of 
cos 8 = — 1.19. These two values are obtained for the best fit values of sin 2 023 and sin 2 0 i 3 . 
However, for sin 2 023 = 0.419 we find the physical value cosd = —0.990, while in the part of 
the 3 ct allowed interval of sin 2 023, 0.374 < sin 2 023 < 0.419, we have —0.495 > cos 8 > —0.990. 

If Gf = A 3 , we find phenomenologically viable results i) for G e = Z 3 , in the case Bl, 
ii) for G e = Z 3 , in the cases Bl and B2, and iii) for G e = Z 2 x Z 2 , in the case B2. More 
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specifically, if G e = Z 3 , we obtain in the case B1 (sin 2 0^ 2 ; sin 2 0 23 ) = (1/3,1/2) leading to 
the predictions sin 2 $12 = 0.341 and cos 5 = 0.570. For G e = Z 5 in the case B1 (case B2) 
we find (sin 2 0° 2 , sin 2 $ 23 ) = (0.276,1/2) ((sin 2 0° 2 , sin 2 6 ^ 3 ) = (0.138,0.160)), which leads to 
the predictions sin 2 $12 = 0.283 and cos 5 = 0.655 (sin 2 @12 = 0.259 and cos <5 = —0.229). 
Finally, for G e = Z 2 x Z 2 in the case B2 we have two sets of values for (sin 2 0 o y2 . sin 2 # 33 ). The 
first one, (sin 2 0\ 2 i sin 2 0J 3 ) = (0.096,0.276), together with the best fit values of sin 2 #13 and 
sin 2 6 * 23 , leads to sin 2 0 12 = 0.330 and cos <5 = —1.36. However, cos 5 takes the physical value of 
cos 5 = —0.996 for sin 2 023 = 0.518. In the part of the 3a allowed interval of values of sin 2 623 , 
0.518 < sin 2 6*23 £ 0.641, we have —0.996 < cos S < —0.478. For the second set of values, 
(sin 2 ^ 2 , sin 2 ^ 3 ) = (1/4,0.127), we get the predictions sin 2 $12 = 0.330 and cosd = 0.805. 


5 The Pattern G e = Z 2 and G v = Z 2 


In this section we derive sum rules for cos 5 in the case given in eq. (24). We recall that 


when the residual symmetries are G e = Z 2 and G v = Z 2 , each of the matrices U e and U v is 
determined up to a complex rotation in one plane. 


5.1 The Case with U12G12Z' 
Cl) 


and [/i 3 (6*3 3 , £33) Complex Rotations (Case 

Similar to the already considered cases we combine the parametrisation of the PMNS matrix 


U given in eq. (24) with (ij) = (12) and (rs) = (13), with the parametrisation of U° given as 

u°{ei 2 , 0 1 ° 3 , e° 23 , s° 12 , si 3 ) = u 12 (o° 12 , s^R^e^u^, <sy, (7i) 

and get the following expression for U (as usual, we refer to Appendix [B] for details): 

U = U 12 w 2 , St2)U l 2{e° 1 2, S° 12 )R23(e°23)Ui3(9° 13 , S° vs )Ur3(9G %)Qo • (72) 


Utilising the results derived in Appendix [B| and reported in eq. (212), we can recast eq. 
in the form 


U = Ri2(e e 12 )Pi(6)R23(9° 2 3)Ri3(9^)Q 0 , PZS) = diag(e* 5 ,1,1). 


(73) 


Here 5 = a e — /3 e + a u + /3 U and we have redefined the matrix Q o by absorbing the diagonal 
phase matrix Pi3(—/3 l ',—a u ) = diag(e _ *^ I/ , 1, e~ ia '') in it. Using eq. (73) and the standard 
parametrisation of the PMNS matrix U, we finch 

sin 2 6*13 = | C e 31 2 = cos 2 0\ 2 sin 2 0" 3 + cos 2 0\ 3 sin 2 0\ 2 sin 2 6* 23 


+ — sin 20^2 sin 20^ sin 0 2 3 cos <5 , 


sin 2 023 = 
sin 2 0 i 2 = 


I U, 


fj,3\ 


sin 2 0^3 — sin 2 033 + cos 2 6” 3 sin 2 0 2: 


13 ■ 


’23 


1 - |c e3 | 2 

I £421 2 


1 — sin 2 0 i 3 


sin 2 0®2 cos 2 0 23 


(74) 

(75) 

(76) 


1 - |Ce 3| 2 1 -sin 2 013 

The sum rule for cos 5 of interest can be derived by comparing the expressions for the 
absolute value of the element U T 2 of the PMNS matrix in the standard parametrisation and 


in the one obtained using eq. (73): 


I U T 2 1 = | cos 0 i 2 sin 023 + sin 0 i 3 cos 0 2 3 sin 0 i 2 e“ | = | sin 0. 


231 


(77) 
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From the above equation we get for cos 5: 


cos 5 


sin 2 $23 — cos 2 0 i 2 sin 2 023 — cos 2 623 sin 2 0 i 2 sin 2 0 i 3 
sin 0 i 3 sin 2023 sin 0 i 2 cos 0 i 2 


(78) 


Given the assumed breaking pattern, cos S depends on the flavour symmetry Gf via the value 
of 02 3 . Using the best fit values of the standard mixing angles for the NO neutrino mass 
spectrum and the requirement | cos <7| < 1 , we find that sin 2 0^3 should lie in the following 
interval: 0.236 < sin 2 0^3 < 0.377. Fixing two of the three angles to their best fit values and 
varying the third one in its 3<r experimentally allowed range and considering all the three 
possible combinations, we get that | cos <5| < 1 if 0.195 < sin 2 0^3 < 0.504. 


5.2 The Case with Ui 3 (9l 3 ,5f 3 ) an d ^ 12 (^ 12 )^ 12 ) Complex Rotations (Case 
C2) 


As in the preceding case, we use the parametrisation of the PMNS matrix U given in eq. (24) 
but this time with (ij) = (13) and (rs) = (12), and the parametrisation of U° as 


to get for U (again the details can be found in Appendix |B| : 


The results derived in Appendix [B] and reported in eq. ( 212 ) allow us to rewrite the expression 
for U in eq. (80) as follows: 


^ 12 )) 

(79) 

S u 12 )Qo. 

(80) 


U = Ri 3 ( 9 e 13 )P 1 ( 8 )R 23 ( 6 o 23 )R 12 (e^)Q 0 , P^S) = diag(e i5 , 1,1), 


(81) 


where 5 = a e — f3 e + oP + f3 u , and also in this case we have redefined the matrix Q$ by 
absorbing the phase matrix Pi 2 (—( 3 u ,—a L ') = diag(e - *^, e~ ia ’', 1) in it. From eq. (81) and 
the standard parametrisation of the PMNS matrix U we get: 

sin 2 0 i 3 = | U e 3 1 2 = cos 2 023 sin 2 0 ^ 3 , 


sin 2 023 = 




sin 2 0 ^ 


sirr 012 = 


l-|£7 e3 | 2 cos 2 0 i 3 ’ 
|C e2 | 2 1 


1 | C e 31 2 1 -sin 2 0 i 3 


[ cos 2 0 i 3 sin 2 0 i 2 + cos 2 0 i 2 sin 2 0 f 3 sin 2 0 ^ 


— - sin 2013 s i n 2^12 s i n ^23 cos <5] . 


(82) 

(83) 

(84) 


Given the value of sin 2 0 ^ 3 , eq. (83) implies the existence of a correlation between the values 
of sin 2 023 and sin 2 0 i 3 . 

Comparing the expressions for the absolute value of the element LA 1 of the PMNS matrix 
in the standard parametrisation and in the one obtained using eq. (81), we have 


| Up 11 = | sin 0 i 2 cos 023 + sin 0 i 3 sin 0 23 cos 0 i 2 e* 5 1 = | sin 0 32 cos 2 0^ 3 1 . 


(85) 


19 













From the above equations we get for cos 5: 


- cos 2 013 (cos 
coso =- 


,2 /jo 


In this case cos 5 is a function of the known neutrino mixing angles 6 12 and @ 13 , of the angle 0 23 
fixed by Gf and the assumed symmetry breaking pattern, as well as of the phase parameter 
5 of the scheme. Predictions for cos 5 can only be obtained when 5 is fixed by additional 
considerations of, e.g., GCP invariance, symmetries, etc. In view of this we show in Fig. [l] 
cos 8 as a function of cos 8 for the current best fit values of sin 2 9 12 and sin 2 $ 13 , and for the 
value sin 2 023 = 1/2 corresponding to Gf = £ 4 . We do not find phenomenologically viable 
cases for A 4 (T') and A; 5 . Therefore we do not present such a plot for these groups. 

cos 8 



Figure 1: Dependence of cosd on cosd in the case of Gf = S 4 with sin 2 $23 = 1/2- The 
mixing parameters sin 2 $12 and sin 2 $13 have been fixed to their best fit values for the NO 
neutrino mass spectrum quoted in eqs. © and Q. The solid (dashed) line is for the case 
when sin 2 0f 3 sin 26 1 2 is positive (negative). 


5.3 The Case with ^12(^12, ^1^) an d U 23 (9 23 ,S 23 ) Complex Rotations (Case 
C 3 ) 

We get for the PMNS matrix U, 

u = U 12 (0t 2 , 5 e 12 )U 12 (9° 12 ,8° 12 )R 13 (9° 13 )U 23 (6° 23 ,8° 23 )U 23 (9" 23 ,5" 23 )Q 0 , (87) 

utilising the parametrisations of U shown in eq. ( |24[ ) with (ij) = (12) and (rs) = (23) and 
that of U° given below (further details can be found in Appendix [bJ) , 

U°(e° 12 ,e° 13 ,0° 23 ,8° 12 ,5° 23 ) = Uu(01 2 ,51 2 )R 13 (91 3 )U 23 (9° 23 ,5° 23 ) . ( 88 ) 
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With the help of the results derived in Appendix [B| and especially of eq. (212), the expression 
in eq. (87) for the PMNS matrix U can be brought to the form 


U = Ri 20 t 2 )P 2 ( 6 )R 13 ( 6 ° 13 )R 23 (O" 23 )Q o , P 2 (S) = diag(l, e iS , 1), (89) 


where 6 = (3 e — a e + a v + (5 V and, as in the preceding cases, we have redefined the phase 
matrix Qo by absorbing the phase matrix P 2 3 (—P 1 ', —a u ) = diag(l, e~ l ^ v , e~ iaV ) in it. Using 
eq. (89) and the standard parametrisation of the PMNS matrix U, we find: 


sin 2 0 i 3 = | C/ e 3 1 2 = sin 2 0\ 2 sin 2 0 23 + cos 2 0\ 2 cos 2 8 23 sin 2 9l 3 

+ — sin 2012 s i n ^9 23 sin 0i 3 cos <5 , 

2 Q I Ufj.3 1 2 sin 2 023 — sin 2 0 i 3 + cos 2 0o 3 sin 2 d\ 3 

5111 = rVKSP = 1 - sin 2 0i3 1 

■ 2/1 I U e 2 1 2 sin 2 0f 2 — sin 2 0 i 3 + cos 2 0? 2 sin 2 0? 3 

sm 012 = 1-17712 = -I--- 

1 - \Ue3\ 1 — sm 2 013 


(90) 

(91) 

(92) 


The sum rule for cos 5 of interest can be derived, e.g., by comparing the expressions for the 
absolute value of the element U T i of the PMNS matrix in the standard parametrisation and 
in the one obtained using eq. (89): 


\U T \\ = | sin0i2 sin 9 23 — sin0i3 cos0i2 cos d 23 e l5 1 = | sin0° 3 


(93) 


For cos 5 we get: 


cos 6 = 


sin 2 0i2 shr 023 — sin 2 0 33 + cos 2 0i2 cos 2 9 23 sin 2 0 i 3 
sin 0 i 3 sin 2023 sin 0 i 2 cos 0 i 2 


(94) 


In this case, in contrast to that considered in the preceding subsection, cos 6 is predicted 
once the angle 0J 3 , i.e., the flavour symmetry Gf, is fixed. Using the best fit values of 
sin 2 0i2, sin 2 0i3 and sin 2 023 for the NO neutrino mass spectrum, we find that physical values 
of cos 6 satisfying | cos 5| < 1 can be obtained only if sin 2 0J 3 lies in the following interval: 
0.074 < sin 2 0^ < 0.214. Fixing two of the three neutrino mixing parameters sin 2 0i2, sin 2 0i3 
and sin 2 9 23 to their best fit values and varying the third one in its 3cr experimentally allowed 
range and taking into account all the three possible combinations, we get that |cos<5| < 1 
provided 0.056 < sin 2 0 33 < 0.267. 


5.4 The Case with ^ 13 ( 6 * 13 , 5i 3 ) and U 23 (9 23 ,5 23 ) Complex Rotations (Case 

Cl) 

The parametrisation of the PMNS matrix U, to be used further, 

U = Ci 3 (0?3, 6 e i3)Ui3(9° 13 ,61^2(012^23(8^, 6 ^ 23 ( 8 ^ 3 , SZ 3 )Qo , (95) 

is found in this case from the parametrisations of the matrix U given in eq. ( |24[ ) with (ij) = 
(13) and (rs) = (23) and that of U° shown below (see Appendix [B| for details), 

U °( 0 i 2 , 9 ° l3 , 023, 5°3, (^23) = Ui3(0i3, ^13)^12(^12)^23(^23’ ^23) ■ (^6) 
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The results presented in eq. (212) of Appendix |B| allow us to recast eq. ( f95| in the form: 

U = R 13 (e e 13 )P 3 (S)R 12 (e o 12 )R 23 (^ 3 )Q 0 , P 3 (5) = diag(l, 1, e"). (97) 

Here 5 = (3 e —a e —a u —(3 l/ and we have absorbed the phase matrix P 23 (a 1 ', f3 u ) = diag(l, e taV , e 1 ^’ 
in the matrix Q q. Using eq. (97) and the standard parametrisation of the PMNS matrix U, 
we find: 


sin 2 % 


sin 2 023 
sin 2 0 i 2 


= | t/ e 31 2 = cos 2 023 sin 2 § 1 3 + cos 2 0 i 3 sin 2 023 sin 2 0\ 2 

H— sin 26 e l3 sin 20^ sin 0^ 2 cos <5 , 

I Ufi 3 1 2 cos 2 0J 2 sin 2 023 
1 - \U e3 \ 2 1-sin 2 0i3 

| f/ e 2 1 2 _ sin 2 §l 3 — sin 2 0 i 3 + cos 2 0f 3 sin 2 0J 2 

1 - \Ue 3 \ 2 1-sin 2 0 i 3 


(98) 

(99) 

( 100 ) 


Comparing the expressions for the absolute value of the element U^i of the PMNS matrix in 
the standard parametrisation and in the one obtained using eq. (97), we find 


\Ufj,i \ = | sin 0 i 2 cos 023 + sin 0 i 3 sin 023 cos 0 i 2 e 115 1 = | sin0J 2 1 
From the above equation we get for cos 5: 

sin 2 0\ 2 — cos 2 023 sin 2 0 i 2 — cos 2 0 i 2 sin 2 0 i 3 sin 2 0 23 


cos 5 = 


sin 0 i 3 sin 2023 sin 0 i 2 cos 0 i 2 


( 101 ) 


( 102 ) 


The predicted value of cos 6 depends on the discrete symmetry G / through the value of the 
angle 0J 2 . Using the best fit values of the standard mixing angles for the NO neutrino mass 
spectrum and the requirement | cos <7| < 1, we find that sin 2 0^ 2 should lie in the following 
interval: 0.110 < sin 2 0^ 2 < 0.251. Fixing two of the three neutrino mixing angles to their best 
fit values and varying the third one in its 3cr experimentally allowed range and accounting for 
all the three possible combinations, we get that | cos<5| < 1 if 0.057 < sin 2 0\ 2 < 0.281. 


5.5 The Case with U 23 (9 23 ,5 23 ) and Ui 3 (9^ 3 , 5^) Complex Rotations (Case 
C5) 

The parametrisation of the PMNS matrix U, which is convenient for our further analysis, 


u = U 23 (d 2 e 3, S e 23 )U 23 (0° 23 ,6 23 )R\2 (0°2) U± 3 (0i3, <5 i ° 3 )17 i 3 (0?3, $rs)Qo 


(103) 


can be obtained in this case utilising the parametrisations of the matrix U given in eq. ( |24| ) 
with (ij ) = (23) and (rs) = (13) and that of the matrix U° given below (for details see 
Appendix [B|), 


U°(9l 2, 0i3 1 02°3, <5i°3, ^ 23 ) = U 23 (0° 23 ,6° 23 )R 12 (0° 12 )U 13 (e° 13 , S° 13 ) 


(104) 


The expression in eq. 
eq. (212) in Appendix |Bl 


for U can further be cast in a “minimal” form with the help of 


u = R 23 (e e 23 )P 3 (5)R 12 (e° 12 )R 13 (e u 13 )Q 0 , p 3 (5) = diag(l, 1, e"), (105) 
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where 5 = (3 e — a e — a v — /3 U and we have absorbed the matrix Pi 3 (a", /?") = diag(e“^, 1, e 1 ^) 
in the phase matrix Qq. Using eq. (105) and the standard parametrisation of the PMNS matrix 
U, we find: 


sin 2 


\U „ 3 


n _ I I _ _ 

23 1 | U e 31 2 1 — sin' 

1 


2 013 


[ cos 2 013 sin 2 023 




- sin 20 


■2 a I Ue2 | 2 sill 2 9° l2 

sm 9,2 = r^F = i- 


*23 sin 2013 s i R ^12 cos ) 


— sin 


2 013 


(106) 


(107) 

(108) 


We note that, given Gf, the values of sin 2 0 i 2 and sin 2 0 i 3 are correlated. This allows one to 
perform a critical test of the scheme under study once the discrete symmetry group Gf has 
been specified. 

The sum rule for cos <5 of interest can be derived, e.g., by comparing the expressions for 
the absolute value of the element U T 2 of the PMNS matrix in the standard parametrisation 
and in the one obtained using eq. (105): 


\U t2 \ = | cos 0 i 2 sin 023 + sin 0 i 3 cos 023 sin 0 i 2 e* 5 j = | cos 0 ° 2 sin 1 


(109) 


This leads to 

cos 2 0 i 3 (cos 2 0 i 2 sin 2 0 f 3 — sin 2 @ 23 ) + sin 2 0 j 2 (sin 2 023 — cos 2 623 sin 2 0 i 3 ) 


cos 8 = 


sin 2023 sin 0 i 3 | sin 0 J 2 1 (cos 2 0 i 3 — sin 2 0 ° 2 ) 2 


( 110 ) 


Similar to the case C2 analysed in subsection [572j cos <5 is a function of the known neutrino 
mixing angles 0 i 3 and 023 , of the angle 0\ 2 fixed by Gf and the assumed symmetry breaking 
pattern, as well as of the phase parameter 5 of the scheme. Predictions for cos 6 can be 
obtained if 6 is fixed by additional considerations of, e.g., GCP invariance, symmetries, etc. 
In view of this we show in Fig. [ 2 ] cos <5 as a function of cos <5 for the current best fit values of 
sin 2 0 i 3 and sin 2 023, and for the value sin 2 0 \ 2 = 1/4 corresponding to Gf = S 4 and A$. We 
do not find phenomenologically viable cases for A 4 ( T '). Therefore we do not present such a 
plot for these groups. 


5.6 The Case with U 23 (9 23 , S 23 ) and ^ 12 (^ 12 )^ 12 ) Complex Rotations (Case 
C6) 

We show below that in this case cos 5 coincides (up to a sign) with the cosine of an uncon¬ 
strained CPV phase parameter of the scheme and therefore cannot be determined from the 
values of the neutrino mixing angles and of the angles determined by the residual symme¬ 
tries. Indeed, using the parametrisation of the matrix U given in eq. ( |24| ) with (ij ) = (23) 
and (rs) = (12) and the parametrisation of U° as follows (see Appendix |B| for details), 

U°(9° l2 ,0° 13 ,9° 23 ,5° 12 ,5° 23 ) = U 23 (0° 23 ,6° 23 )R 13 (9l 3 )U 12 (01 2 ,5° 12 ), (111) 

we get for U: 

U = C 23 (0| 3 , 5 e 23 )U 23 (9° 23 ,6° 23 )R 13 (9° 13 )U 12 (0? 2 , S° 12 )U 12 (^ 2 , % 2 )Q 0 . ( 112 ) 
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cos <5 



Figure 2: Dependence of cosd on cos 5 in the case of Gf = S 4 or A 5 with sin 2 9 32 = 1/4. 
The mixing parameters sin 2 023 and sin 2 $13 have been fixed to their best fit values for the NO 
neutrino mass spectrum quoted in eqs. ([7]) and The solid (dashed) line is for the case 
when sin 20 23 sin 20 33 is positive (negative). 


The results derived in Appendix [B] in eq. (212) make it possible to recast eq. (112) in the 
form: 

U = R2 3 (e e 23 )P 2 (5)Ri 3 (0° 13 )R 12 (0 , { 2 )Qo , P 2 (S) = diag(l, e" 1). (113) 


Here 5 = a e — (3 e — a u — /3 U and, as in the preceding cases, we have redefined the phase matrix 
Qo by absorbing the phase matrix Pi- 2 (a u , fi v ) = diag(e' 
the standard parametrisation of the PMNS matrix U, we find: 


ia \e 1 ^ 1 1) in it. Using eq. (113) and 


sin 2 013 = 
sin 2 023 = 

sin 2 0 i 2 = 


|t / e3 | 2 = sin 2 0 ? 3 , 

I Ufjft | 2 . 2 ne 

= sm 0 23 , 


1-I^e3| 2 

i-|u e3 | 2 


= sin 2 0 ,J 


12 • 


(114) 

(115) 

(116) 


Comparing the absolute value of the element U T \ allows us to find that cos <5 = ±cos<5. 

It follows from eq. (114) that for a given flavour symmetry G/, the value of sin 2 0i 3 is 
predicted. This allows to test the phenomenological viability of the case under discussion, 
since the value of sin 2 0 i 3 is known experimentally with a relatively high precision. 


A comment, analogous to those made in similar cases considered in subsections 3.3 and 
4.3, is in order. Namely, for a non-Abelian flavour symmetry Gf which allows to reproduce 
correctly the observed values of sin 2 0 i 2 , sin 2 0 i 3 and sin 2 9 23 , it might be possible to obtain 
physically viable prediction for cos 6 by employing GCP invariance in the charged lepton or 
the neutrino sector. 
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5.7 The Case with Ui 2 (9l 2 ,8f 2 ) and U\ 2 [6 32 ,8 32 ) Complex Rotations (Case 
C7) 

Using the following parametrisation of U°, 


U°(0° 12 ,e° 12 ,8° 23 ,5 o 12 ,8 o 12 ) = U 12 (e° 12 , 8 ° 12 )R 23 (d° 23 )U 12 (6 

we have for U: 

U = C 1 2(^2^ 1 e 2)t/ 1 2W 2) ^2)^23(02 O 3)^12W2,^2)^12(0 

Utilising the results derived in Appendix |B| and reported in eq. ( |212 ) 
in the form: 


^12)) 

(117) 

% 2 )Qo- 

(118) 


, we can recast eq. ( 118|) 


U = i?i2(e? 2 )Ti(<5)^2 3 (02 O 3)^i2(^2)Qo , Pi(S) = diag(e" 1,1). (119) 


Here 8 = a e — f3 e + a u + f3 y and we have redefined the matrix Q 0 by absorbing the diagonal 
phase matrix Pi 2 (—/3 u ,—a u ) = diag(e _ *^‘",e -1 " 1 ", 1) in it. Using eq. (119) and the standard 
parametrisation of the PMNS matrix U, we finch 

sin 2 0i3 = |U e3 | 2 = sin 2 0 23 sin 2 0{ 2 , 

1 U /t 31 2 sin 2 $23 cos 2 0 f { 2 


sin 2 023 = 


sin 2 0i2 = 


1 — sin 2 0i3 
1 


1 - |C e3 | 2 

We2? _ 

1 - |C e3 | 2 1-sin 2 013 

1 


[ cos 2 0 23 cos 2 0i 2 sin 2 §f 2 + cos 2 0f 2 sin 2 0/ 2 


+ - sin 20 i 2 sin 20 32 cos 0 23 cos 


S\ 


From eqs. (120) and (121) we see that the angles 0i 3 and d 23 are correlated: 


sin 2 0 2 3 — sin z 0i 3 


sin 2 023 = 


( 120 ) 

( 121 ) 

( 122 ) 

(123) 


1 — sin 2 0i3 

Comparing the expressions for the absolute value of the element U T \ of the PMNS matrix 


in the standard parametrisation and in the one obtained using eq. (119), we have 
| U T i | = | sin 0i2 sin 0 23 — sin di 3 cos 0 23 cos di 2 e lS \ = | sin 0^ 2 sin 0 23 1. 
From the above equations we get for cos 5: 


cos 5 = 


sim 0i3(cos 2 0i2 cos 2 0 23 — sin 2 6 \ 2 ) + sin 2 0 23 (sin 2 0i2 — cos 2 0i 3 sin 2 0i 2 ) 
sin 2012 sin 0i31 cos 02 3 |(sin 2 0 23 — sin 2 013)2 


(124) 


(125) 


In this case cos 5 is a function of the known neutrino mixing angles 0i2 and 0i 3 , of the angle 0 23 
fixed by Gf and the assumed symmetry breaking pattern, as well as of the phase parameter 
8 of the scheme. Predictions for cos 8 can only be obtained when 8 is fixed by additional 
considerations of, e.g., GCP invariance, symmetries, etc. In view of this we show in Fig. [3] 
cos 8 as a function of cos 8 for the current best fit values of sin 2 0i2 and sin 2 0i 3 , and for the 
value sin 2 0 23 = 1/2 corresponding to Gf = S4. We do not find phenomenologically viable 
cases for Gf = A 4 (T') and A$. 
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cos <5 



Figure 3: Dependence of cosd on cos<5 in the case of Gf = 64 with sin 2 # 23 = 1/2- The 
mixing parameters sin 2 $12 and sin 2 $13 have been fixed to their best fit values for the NO 
neutrino mass spectrum quoted in eqs. ([d]) and The solid (dashed) line is for the case 
when sin 2 0 f 2 sin 2 #i 2 is positive (negative). 


5.8 The Case with Ui 3 (9f 3 , 5 ® 3 ) and ^ 13 (^ 13 ,^ 13 ) Complex Rotations (Case 
C8) 

Using the following parametrisation of U°, 

u°{ 6 l 3 , e\ 3 , # 2 ° 3 , s° 13 , 5 °, 3 ) = um 3 , 4 ), ( 126 ) 

we have for U: 

U = UM 3 , 6 e 13 )U 13 (6° 13 ,5 0 V3 )R 23 (e? 23 )U V3 (6 0 V3 ,6° 13 )U 13 ^ 3 , 5" 13 )Q 0 . (127) 

Utilising the results derived in Appendix |B| and reported in eq. ( |212 ), we can recast eq. ( 127| ) 
in the form: 


U = 7?i3(^ 3 )A(5) J R23(02 O 3)^i3(^3)Qo , Pi(S) = diag(e i 5 ,1, 1) . 


(128) 


Here <5 = a e — /3 e + a v + f3 u and we have redefined the matrix Q o by absorbing the diagonal 
phase matrix P\ 3 {—j3 v ,—Q 1 ') = diag(e~*^, 1, e~ ia '') in it. Using eq. (128) and the standard 
parametrisation of the PMNS matrix U, we find: 

Sill 2 #13 = | U e 3 1 2 = COS 2 #23 C °S 2 #13 s bl 2 #13 + COS 2 #13 Sill 2 #43 


sin 2 #23 = 


sin 2 #12 = 


+ - sin 2 #i 3 sin 2 #i 3 cos #23 cos 5 , 
1 ^ 31 2 _ sin 2 #23 cos 2 #13 


i-|c e3 | 2 

e2 


| U e 2 1 2 sin 2 #? Q sin 2 #? 


1 — sin 2 #13 

23 sin2 #13 


i-|c e3 | 2 


1 — sin 2 #13 


(129) 

(130) 

(131) 
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The sum rule for cos 5 of interest can be derived by comparing the expressions for the 
absolute value of the element 2 of the PMNS matrix in the standard parametrisation and 
in the one obtained using eq. (128): 


t/ /i2 1 = | cos 612 cos 023 — sin 0 i 3 sin 023 sin | = | cos 0 23 | 
From the above equation we get for cos 5: 


cos 5 = 


cos 2 0 i 2 cos 2 023 — cos 2 0 23 + sin 2 0 i 2 sin 2 623 sin 2 0 i 3 
sin 0 i 3 sin 2023 sin 0 i 2 cos 0 i 2 


(132) 


(133) 


Given the assumed breaking pattern, cos S depends on the flavour symmetry Gf via the value 
of 0^. Using the best fit values of the standard mixing angles for the NO neutrino mass 
spectrum and the requirement |cos5| < 1, we find that sin 2 0 23 should lie in the following 
interval: 0.537 < sin 2 0 23 < 0.677. Fixing two of the three angles to their best fit values and 
varying the third one in its 3<r experimentally allowed range and considering all the three 
possible combinations, we get that | cos<5| < 1 if 0.496 < sin 2 0 23 < 0.805. 


5.9 The Case with U 2 3(0 23 ,5 23 ) an d ^ 23 (^ 93 , <^ 23 ) Complex Rotations (Case 
C9) 

Using the following parametrisation of U°, 


u°(0°2 3 , 0 2 ° 3 , ei 2 , <5 2 ° 3 , ^ 23 ) = «0 2 O 3 , s° 23 )r 12 (9 0 12 )u 23 (9 0 23 , s° 23 ) , 

we have for U: 

u = U 23 (9 e 23 , ^3)^23(0235 5 o 23 )R 12 (0° 12 )U 23 (e o 23 , S? 23 )U 23 (9^ 5" 23 )Q 0 


(134) 


(135) 


Utilising the results derived in Appendix |B| and reported in eq. ( |212 ), we can recast eq. ( 135| ) 
in the form: 


u = R23(9 e 23 )P2(d)R 1 2(9 o 12 )R23(9 u 23 )Q 0 , P 2 (S) = diag(l, e ia , 1). 


(136) 


Here 6 = a e — (3 e + a u + and we have redefined the matrix Qq by absorbing the diag¬ 
onal phase matrix P 2 3 {a u 1 l 5 u ) = diag(l, e ia '', e 1 ^) in it. Using eq. (136) and the standard 
parametrisation of the PMNS matrix U, we find: 


sin 2 0 i 3 = |U e3 | 2 = sin 2 0i 2 sin 2 0 23 , 

I 1 2 _ 1 


sin 0 23 = 


1 - |C e3 | 2 1-sin 2 0 i 3 


[ cos 2 0i 2 cos 2 0| 3 sin 2 0 23 + cos 2 0 23 sin 2 9 23 


1 


+ - sin 2 9 23 sin 20 23 cos 0J 2 cos 


S\ 


sin 2 0i 2 = 


|U e2 | 2 


sin 2 91 2 cos 2 0 23 


1 - |C e3 | 2 1-sin 2 0 i 3 

From eqs. (137) and ( |139 ) we find that the angles 0i3 and 0i2 are correlated: 

sin 2 91 2 — sin 2 0 i 3 


siir 0i2 = 


1 - sin 2 0 


(137) 

(138) 

(139) 

(140) 
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Comparing the expressions for the absolute value of the element U T \ of the PMNS matrix 
in the standard parametrisation and in the one obtained using eq. (136), we have 


\U T i\ = | sin0i2 sin 023 — sin0i3cos023cos0i2e l<5 | = | sin0 23 sin0 32 1 • 

From the above equations we get for cos 5: 

sin 2 0 i 3 (cos 2 023 cos 2 0i 2 — sin 2 023 ) + sin 2 0 32 (sin 2 023 — cos 2 0 i 3 sin 2 0| 3 ) 
sin 2023 sin 0 i 3 1 cos 0J 2 |(sin 2 0 32 — sin 2 013 ) 2 

In this case cos 5 is a function of the known neutrino mixing angles 023 and 0i3, of the angle 0^ 2 
fixed by Gf and the assumed symmetry breaking pattern, as well as of the phase parameter 
5 of the scheme. Predictions for cos <5 can only be obtained when 5 is fixed by additional 
considerations of, e.g., GCP invariance, symmetries, etc. In view of this we show in Fig. [4] 
cos 5 as a function of cos 5 for the current best fit values of sin 2 023 and sin 2 0i3, and for 
the value sin 2 0 32 = (r + 2)/(4r + 4) = 0.345 corresponding to Gf = A 5 . We do not find 
phenomenologically viable cases for Gf = A4 ( T') and S4. 


(141) 

(142) 


cos 8 



Figure 4: Dependence of cos 5 on cos 8 in the case of Gf = A$ with sin 2 0 32 = (r+2)/(4r+4) = 
0.345. The mixing parameters sin 2 0 23 and sin 2 0i3 have been fixed to their best fit values for 
the NO neutrino mass spectrum quoted in eqs. Q and Q. The solid (dashed) line is for the 
case when sin 20f 3 sin 20 23 is positive (negative). 


5.10 Results in the Cases of Gf = A 4 ( T '), S 4 and A 5 


The schemes considered in Sections 5.1 — 15.9| can be applied when considering the breaking 
Gf to G e = Z2 and G v = Z2, for both Majorana and Dirac neutrinos. As explicit examples 
of this, we now consider Gf = A4 ( T '), S4 and A 5 broken to G e = Z2 and G v = Z2. As such, 
we have considered all possible combinations of residual Z 2 symmetries for a given flavour 
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symmetry group, namely, G e = Z 2 and G v = Z 2 for Gf = A4 ( T'), S4, A5. For instance, in 
the cases of the schemes described in subsections 5.1 - |5.5[ and Gf = S4 broken to G e = Z% 
and G v = Z\ with (a, 6) = ( T 2 U,U ), ( T 2 U,SU ), (T^U\TU), ( T 2 U,STSU ), etc. (a total 
of 24 combinations of order two elements), the value of the relevant parameter contained in 
the fixed matrix U° yields sin 2 6^3 = 1/4, sin 2 ^3 = 1/2, sin 2 = 1/4, sin 2 0 ° 2 = 1/4, and 
sin 2 OI 2 = 1/4, respectively. In A§ for the cases Cl, C3, C4 and C5 we find the sine square 
of the corresponding fixed angle in the matrix U° to be 1/4, e.g., for G e = Z% and G u = Z\ 
with (a, b) = (S, ST 2 ST 3 S), ( S,ST 3 ST 2 S ), ( S,T 2 ST 3 ), ( S,T 3 ST 2 ), etc. (in total, for 60 
combinations of order two elements). 

For the symmetry group A4 we find that none of the combinations of the residual sym¬ 
metries G e = Z 2 and G u = Z 2 provide physical values of cos 5 and phenomenologically viable 
results for the neutrino mixing angles simultaneously. 

For Gf = S 4 , using the best fit values of the mixing angles 0 12 , 0 13 and # 23 , we get cos <5 = 
—0.806, —1.52 and 0.992 in the cases Cl, C3 and C4, respectively. Physically acceptable 
value of cos 6 in the case C3 can be obtained for sin 2 023 = 0.562 allowed at 3a, for which 
cos 5 = —0.996. In the part of the 3a allowed range of sin 2 023 , 0.562 < sin 2 023 < 0.641, 
we have —0.996 < cosd < —0.690. Further, in the case C2, in which the relevant parameter 
sin 2 0^3 = 1 / 2 , the value of cos 5 is not fixed, while the atmospheric angle is predicted to 
have a value corresponding to sin 2 023 = 0.5 1 2 . Similarly, in the case C5 the value of cos 5 
is not fixed, while sin 2 0 i 2 = 0.256 (which is slightly outside the corresponding 3 a interval). 
In the case C7 we find that cos 5 is not fixed and sin 2 023 = 0.488. Finally, for C 8 with 
sin 2 0^3 = 1/2 and 3/4, using the best fit values of the neutrino mixing angles for the NO 
spectrum, we have cos <5 = —1.53 and 2.04, respectively. The physical values of cos 8 can be 
obtained, using, e.g., the values of sin 2 023 = 0.380 and 0.543, for which cos 5 = —0.995 and 
0.997, respectively. In the parts of the 3a allowed range of sin 2 823 , 0.374 < sin 2 623 < 0.380 
and 0.543 < sin 2 023 < 0.641, we have —0.938 > cos <5 > —0.995 and 0.997 > cos <5 > 0.045, 
respectively. 

For the A$ symmetry group the cases Cl with sin 2 0o 3 = 1/4, C3 with sin 2 0 J 3 = 1/4 
and C4 with sin 2 0 ^ 2 = 1/4 lead to the same predictions obtained with Gf = S4, namely, 
cos 6 = —0.806, —1.52 and 0.992, respectively. Moreover, in the case C3 (case C4) the 
value of sin 2 0^ 3 = 0.096 (sin 2 0 ° 2 = 0 .096) is found, which along with the best fit values of 
the mixing angles gives cos 5 = 0.688 (cos<5 = —1.21). Using the value of sin 2 023 = 0.487 
allowed at 2a, one gets in the case C4 cosd = —0.997, while in the part of the 3a allowed 
range of sin 2 023 , 0.487 < sin 2 023 < 0.641, we have —0.997 < cos 6 < —0.376. Note also, if 
sin 2 023 is fixed to its best fit value, one can obtain the physical value of cos 5 = —0.999 using 
sin 2 0 i 2 = 0.277. For the part of the 3a allowed range of sin 2 0 i 2 , 0.259 < sin 2 0 i 2 < 0.277, 
one gets —0.871 > cos 5 > —0.999. The cases C5 and C 8 are the same as for the S 4 symmetry 
group. Finally, in the case C9 the value of cos 6 is not fixed, while using the best fit value of 
the reactor angle, we get sin 2 0 i 2 = 0.330. 


6 Summary of the Results of Sections 3, 4 and 5 

The sum rules derived in Sections 3, 4 and 5 are summarised in Tables [3] and |4j The formulae 
for sin 2 0 i 2 , sin 2 0 i 3 and sin 2 623 , which lead to predictions for the indicated neutrino mixing 
parameters once the discrete flavour symmetry Gf is fixed, are given in Tables [5] and [ 6 j In 
the cases in Tables [5] and [ 6 ] in which cos <5 is unconstrained, a relatively precise measurement 
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of sin 2 0 i 2 , sin 2 $13 or sin 2 023 can provide a critical test of the corresponding schemes due to 
constraints satisfied by the indicated neutrino mixing parameters. 

A general comment on the results derived in Sections 3, 4 and 5 is in order. Since we do 
not have any information on the mass matrices, we have the freedom to permute the columns 
of the matrices U e and U u , or equivalently, the columns and the rows of the PMNS matrix U. 
The results in Tables [3] and [4] cover all the possibilities because, as we demonstrate below, the 
permutations bring one of the considered cases into another considered case. For example, 
consider the case of U = fTi 3 (0f 3 , 5h)u°u 23 (e" 23 , 5% 3 )Qq. The permutation of the second and 
the third rows of U is given by 7 ^ 3 U = '^23 Ui3 {9 33 ,5\ 3 )'K23'^23U°U23(0 2 3-,5 23 )Qo, where we 
have defined 


T23 



(143) 


Since the combination 7 r 23 ^ 13 ( 013 , <5f 3 )7T23 gives a unitary matrix £/i2(0f 3 , <5f 3 ), the result after 
the redefinition, 0f 3 —> 0f 2 , 5f 3 —> b\ 2 and 7723 U° —> U°, yields 


u = u l 2 {ei 2 , 5 t 2 )u°U 23 {e^ 5y 23 )q {) , 

which represents another case present in Table |4j It is worth noting that the freedom in 
redefining the matrix U° follows from the fact that U° is a general 3x3 unitary matrix 
and hence can be parametrised as described in Section [2] and in Appendix [Bj All the other 
permutations should be treated in the same way and lead to similar results. 
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Case 

Parametrisation of U 

Sum rule for cos 5 

Al 

Ul2(9t 2 , St 2 )U 12 (6l 2 , 6° 12 )R 23 (9° 3 )R l3 (9° 13 )Q 0 

cos 2 0i3(sin 2 0 23 — cos 2 0i 2 ) + cos 2 0\ 3 cos 2 0 23 (cos 2 0i2 — sin 2 0i2 sin 2 013) 

sin 2012 sin 0i31 cos0 33 cos0 23 (cos 2 0i3 — cos 2 0 33 cos 2 023)2 

A 2 

Ui 3 (e e l3 , Sh)Ui3(e° 13 , 6° l3 )R 23 (9° 23 )R l2 (9° 12 )Qo 

cos 2 0i3(cos 2 0J 2 cos 2 0 23 — cos 2 0i2) + sin 2 0 23 (cos 2 0i2 — sin 2 0i2 sin 2 013) 

sin 2012 sin 0i 3 sin0 23 (cos 2 0i3 — sin 2 0 23 )2 

A 3 

U 23 ( 0 e 23 , S e 23 )U 23 ( 0 ° 23 , 5° 23 )R 13 (6° 13 )R 12 (e° 12 )Q 0 

± cos 5 2 3 


B1 


^23(^13)^12(^ 2 )^13(^3, ^13)^13(^135 %)Q0 


sin 2023 sin 0i31 sin 0^ 2 1 (cos 2 0i3 


cos 2 0i3(sin 2 0 32 — cos 2 @23) + cos 2 0° 2 cos 2 0J 3 (cos 2 023 — sin 2 0i3 sin 2 023) 
0131 COS 0 32 COS 0131 (COS 2 013 — COS 2 0J 2 COS 2 0J 3 ) 2 


B2 


Rl 3 (0° 13 )R 12 (9° 12 )U 23 (6° 23 , S° 23 )U 23 (9^, <%)Qo 


sin 2023 sin 0i31 


^23^2 O 3)^13W 3 )t / 12(^25^2)^12(^25^2)Q0 ± COS S 12 


B3 


Table 3: Summary of the sum rules for cos <5. The cases Al, A2 and A3 correspond to G e = Z 2 and G u = Z n , n > 2 or Z n x Z m , n, m > 2, 
while Bl, B2 and B3 correspond to G e = Z n , n> 2 or Z n x Z m , n,m > 2 and G v = Z 2 . See text for further details. 




Case Parametrisation of U 


Sum rule for cos 5 


Cl 

Ui 2 (e e u , 

5f 2 )f/ 12 (^ 2 ,^ 2 )i? 23 (0 2 O 3)C/l3W 3 >^3)^13(^3^T3)Qo 

sin 2 0 23 — cos 2 0i 2 sin 2 0 23 — cos 2 0 23 sin 2 0i 2 sin 2 0i 3 


sin 0 i 3 sin 20 23 sin 0i 2 cos 0i 2 

C9 

Ui 3 ( 6 f 3 , 

5f 3 )C/l3W 3 ,%)i?23(0 2 O 3 )^12W 2 ,5? 2 )C/l2(0r2^T 2 )Qo 

cos 2 0i 3 (cos 2 0 23 sin 2 0 32 — sin 2 0i 2 ) + sin 2 0 23 (sin 2 0 32 — cos 2 0i 2 sin 2 0i 3 ) 


sin 20 i 2 sin 0i 3 sin0 23 (cos 2 0i 3 — sin 2 0 23 )2 


Ui 2 ( 0 f 2 , 

6 e i 2 )Ul 2 (e° 12 , S° 12 )Rl3(0 O 13 )U23(0° 23 , 5 2 ° 3 )c/ 23 (^ 3 , ^ 3 )Qo 

sin 2 0i 2 sin 2 0 23 — sin 2 0J 3 + cos 2 0i 2 cos 2 0 23 sin 2 0i 3 


sin 0i 3 sin 20 23 sin 0i 2 cos 0i 2 

Cd 

Ui 3 ( 6 f 3 , 

^13)^13(^1 3 ) ^ 13)^12 (^1 2 )C 2 3 (0 23 )5 23 )C/ 23 (0 23 , S 23 )Qo 

sin 2 0 32 — cos 2 0 23 sin 2 0 i 2 — cos 2 0 i 2 sin 2 0 33 sin 2 0 23 


sin 0 i 3 sin 20 23 sin 0 i 2 cos 0 i 2 


U23(0h, 

^ 23 )^ 23 (^ 23 ) ^ 23 )^ 12 (^ 12 )^ 13 (^ 13 ) ^ 13 )^ 13 (^ 13 ) Wo 

cos 2 0 i 3 (cos 2 0i 2 sin 2 0 23 — sin 2 0 23 ) + sin 2 0 32 (sin 2 0 23 — cos 2 0 23 sin 2 0i 3 ) 


sin 2023 sin 0i 3 sin0 32 (cos 2 0i 3 — sin 2 0 32 )5 

C 6 

U23(0 e 23 , 

^ 23 )^ 23 (^23) ^ 23 )^ 13 (^ 13 )^ 12 (^ 12 ) ^ 12 )^ 12 (^ 12 ) <^12)Qo 

± cos 5 

C7 

Ui 2 (e e 12 , 

5f 2 )f/ 12 (^ 2 ,5^ 2 )i? 23 (0 2 O 3 )^l 2 W 2 ,4)t / 12(^2)^2)Qo 

sin 2 0i 3 (cos 2 0 i 2 cos 2 0 23 — sin 2 0i 2 ) + sin 2 0 23 (sin 2 0i 2 — cos 2 0i 3 sin 2 0 32 ) 


sin 2012 sin 0 i 3 cos0 23 (sin 2 0 23 — sin 2 0 i 3 )2 


Ul3(0f 3 , 

5f3)^13W3)%)^23(02 O 3)^13W3)4)^13(^3)^3)Qo 

cos 2 0 i 2 cos 2 0 23 — cos 2 0 23 + sin 2 0 i 2 sin 2 0 23 sin 2 0i 3 


sin 0 i 3 sin 20 23 sin 0i 2 cos 0i 2 

ro 

U23(0 e 23 , 

^ 23 )^ 23 (^ 23 ’ ^3)^12(^2)^23(023) 4)^23(02^3, 4)^0 

sin 2 0 i 3 (cos 2 0 23 cos 2 0i 2 — sin 2 0 23 ) + sin 2 0 32 (sin 2 0 23 — cos 2 0i 3 sin 2 0 23 ) 


sin 2023 sin 0 i 3 cos0 32 (sin 2 0 32 — sin 2 0 i 3 )2 


Table 4 

: Summary of the sum rules for cos (5. The cases Cl - C9 correspond to G e = Z 2 and G„ = Z 2 . See text for further details. 




Case 

Parametrisation of U 

Sum rule for sin 2 0\ 2 and/or sin 2 0 i 3 and/or sin 2 0 23 

Al 

Ui2(0 e l2 , 5 e n)Un( 0 o 12 ,d o 12 )R 23 ( 8 o 23 )R 13 ( 8 o 13 )Q 0 

sin 2 023 = 

sin 2 0/3 — sin 2 0i 3 + cos 2 0\ 3 sin 2 0 23 

1 — sin 2 0i3 

A2 

^13(0f3^?3)^13(^3^1°3)^23(^3)^12W2)Qo 

sin 2 023 = 

sin 2 0 23 

1 — sin 2 0i3 

A3 

5 e 23 )U 23 ( 8 ° 23 , 6° 23 )R 13 (9° 13 )R 12 (e° 12 )Q 0 

sin 2 0i3 = 

sin 2 0/ 3 , sin 2 0i2 = sin 2 0/ 2 

Bl 

R2 3 (9° 23 )Ri2(e o 12 )Ui3(e° 13 ,61 3 )U 13 W 3 , s^ 3 )q 0 

sin 2 0i2 = 

sin 2 0° 12 

1 — sill 2 013 

B2 

Ri 3 ( 0 ° 13 )Ri 2 ( 0 ° 12 )U 23 (d° 23 , S 23 )U 23 (8 23 , 5 23 )Qo 

sin 2 0i2 = 

COS 2 013 — COS 2 0 32 COS 2 0 ^ 3 

1 — sin 2 0i3 


B3 R 23 (9° 23 )R 13 (e° 13 )U 12 (6° 12 , 6f 2 )U 12 W 2 , S^)Q 0 sin 2 0 13 = sin 2 6° 13 , sin 2 0 23 = sin 2 0 2 ° 3 


Table 5: Summary of the formulae for sin 2 9±2 and/or sin 2 6*13 and/or sin 2 $ 23 - The cases Al, A2 and A3 correspond to G e 
G v = Z n , n > 2 or Z n x Z m , n, m > 2, while Bl, B2 and B3 correspond to G e = Z n , n > 2 or Z n x Z m , n, to > 2 and G„ = Z 2 
for further details. 


= Z 2 and 
. See text 



Case 

Parametrisation of U 

Sum rule for sin 2 9 22 and/or sin 2 @13 and/or sin 2 9 23 

Cl 

Ui 2 (e e 12 , st 2 )Uu(e° 12 ,s° 12 )r 23 (9° 23 )u 13 (9° 13 ,6° l3 )u 13 (9" 13 , 6" 13 )q 0 

not fixed 


C2 

Uis(0i 3 ,6 e 13 )u 13 (e° 13 ,6° 3 )R 23 (o° 23 )u 12 (e° l2 , s o 12 )u 12 (9^,s^ 2 )q 0 

sin 2 023 = 

sin 2 0 ^ 

1 — sin 2 0 i 3 

C3 

Ui 2 (e f 2 , 5l 2 )U 12 (6° 12 , 6° 12 )R 13 (9° 13 )U 23 (9° 23 , 5° 23 )U 23 (9% 3 , S 23 )Qo 

not fixed 


C4 

Ui3(e e 13 , Sf 3 )U 13 (6l 3 ,6° 13 )R 12 (9° 12 )U 23 (9° 23 ,5? 23 )U 23 (6!^ 5" 23 )Q 0 

not fixed 


C5 

U 23 m 3 , Sh)U 23 (e? 23 , 6° 23 )R 12 (0° 12 )Ui 3 (9° 13 , 6l 3 )U 13 ffl 3 , S^Qo 

sin 2 0 i 2 = 

sin 2 0° 2 

1 — sin 2 0 i 3 

C 6 

u 23 (0 e 23 ,8f 23 )U 23 (d? 23 , 5° 23 )Ri 3 (0° 13 )U 12 (0° 12 , 6l 2 )U 12 ffl 2 ,6Z 2 )Qo 

sin 2 0 i 3 = 

sin 2 9 23 

C7 

Ui 2 (9{ 2 , ^i 2 )Ui 2 (9 12 , 5i 2 )R 23 (9 23 )Ui 2 (6i 2 , fii 2 )Ui 2 (9i 2 , Si 2 )Qo 

sin 2 023 = 

sin 2 0^3 — sin 2 0 i 3 

1 — sin 2 0 i 3 

C 8 

Ui 3 (9 e 13 , S e 13 )U 13 (9° 13 ,6° 13 )R 23 (9° 23 )U 13 (9° l3 ,6° 13 )U 13 ffl 3 ,6? 3 )Q 0 

not fixed 


C9 

^23(^23’ ^23)^23(^235 ^23)^12(^12)^23 (9 23 i 5 23 )U 23 (9% 3i S 23 )Qo 

sin 2 0 i 2 = 

sin 2 0 J 2 — sin 2 0 i 3 

1 — sin 2 0 i 3 


Table 6: Summary of the formulae for sin 2 9 12 and/or sin 2 $13 and/or sin 2 $ 23 - The cases Cl - C9 correspond to G e = Z 2 and G v = Z 2 . 
See text for further details. 



7 The Case of Fully Broken G e 


If the discrete flavour symmetry Gf is fully broken in the charged lepton sector the matrix 
U e is unconstrained and includes, in general, three rotation angle and three CPV phase 
parameters. It is impossible to derive predictions for the mixing angles and CPV phases 
in the PMNS matrix in this case. Therefore, we will consider in this section forms of U e 
corresponding to one of the rotation angle parameters being equal to zero. Some of these 
forms of U e correspond to a class of models of neutrino mass generation (see, e.g., m 32 36]) 
and lead, in particular, to sum rules for cos 5. 

We give in Appendix 0 the most general parametrisations of U under the assumption 
that in the case of fully broken G e one rotation angle in the matrix U e vanishes. The second 
case in Table 14 with = 0 have been analysed in 
^ 12 (^ 12 ^ 12 )^ 13 (^ 13 ^ 13 ) h as been investigated in [Ml. 


11 


13 


14 


while the third case with 


7.1 The Scheme with ^23(^23, ^23)^12($12) <^ 2 ) (Case Dl) 

We consider the following parametrisation of the PMNS matrix (see Appendix [Cj first case 
in Table 14): 


U = U 23 (6 e 23 , Sh)Ri2(O 12 )Pi(S)R 23 (0 o 23 )R 13 (0° 13 )Q o , Pi(5) = diag(e , 1,1). (144) 


We find that: 


sin 2 @13 = \U e3 \ 2 = sin 2 6 h 3 ( 0 i 2 ,< 5 , 6 » 33 , 02 3 ), 

IM 2 


sin 2 6*23 = 


sin 2 0 i 2 = 


— sin^ 023(^12) 5 , 023, 623 , 0 33 , 02 3 ) , 


l-|C e3 | 2 

| t/ e 21 2 COS 2 02 3 sin 2 0i2 


(145) 

(146) 

(147) 


l-\U e3 \ 2 COS 2 013 

As it can be seen from the previous equations and the absolute value of the element 2 , 

| Ufj2 1 = | COS 023 COS 012 cos 023 — e ~"^ 3 s j ri s i n $° 3 | . (148) 

a sum rule for cos 5 might be derived in the case of fixed d 23 ■ In the general case of free d | 3 we 
find that cos 5 is a function of b 2 3 - Since in this case the analytical expression of cos 6 in terms 
of b 23 is rather complicated, we do not present this result here. Note that imposing either 


023 = 0 or 0^3 = 0 is not enough to fix the value of cos <5. As eqs. (145) and (146) suggest, 
in the case of fixed 5 23 there exist multiple solutions for the value of cos 5 for any given value 
of b e 23 . This is demonstrated in Fig. [5j in which we plot cos 6 versus b 23 , assuming that the 
angles 0 33 and 0^ 3 have the values corresponding to the TBM, GRA, GRB and HG symmetry 
forms given in Table [ 2 J The figure is obtained for 032 belonging to the first qua drant. The 
solid lines correspond to 6 = cos _1 (cos <5), where cos6 is the solution of eq. (145), whil e the 
dashed lines correspond to 5 = 2ir — cos^ 1 (cos 6 ). Multiple lines reflect the fact that eq. (146) 


for 0| 3 has several solutions. We note that Fig. [H] remains the same for 032 belonging to the 
third quadrant, while for 032 lying in the second or fourth quadrant the solid and dashed lines 
interchange. For the BM (LC) symmetry form cos <5 has an unphysical value, which indicates 
that the considered scheme with the BM (LC) form of the matrix diagonalising the neutrino 
mass matrix does not provide a good description of the current data on the neutrino mixing 
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Figure 5: Dependence of cos S on <5f 3 in the cases of the TBM, GRA, GRB and HG symmetry 
forms. The mixing parameters sin 2 $ 12 , sin 2 d 2 3 and sin 2 do have been fixed to their best fit 
values for the NO neutrino mass spectrum quoted in eqs. (Ji]) - The angle di 2 is assumed 
to belong to the fir st qu adrant. The solid lines correspond to 5 = cos -1 (cos <5), where cos 8 is 
the solution of eq. (145), while the dashed lines correspond to 5 = — cos” 1 (cos <5). See text 

for further details. 


angles 


12 


14 


Thus, we do not present such a plot in this case. If (5| 3 turns out to be fixed 


(by GCP invariance, symmetries, etc.), then, as can be seen from Fig. [5j cos 5 is predicted to 
take a value from a discrete set. For instance, when d |3 = 0 or it , we have 


cos 6 = {—0.135,0.083} for TBM; 
cos 5 = {-0.317,0.269} for GRA; 
cos 5 = {-0.221,0.170} for GRB; 


(149) 

(150) 

(151) 


14 Note that the scheme under discuss ion corresponds to inverse ordering of the charged lepton corrections, 
i.e., u! = U 2 3(e e 2 3 ,S^)Ui2(ef 2 ,Sh) (see ill). 
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cos5 = {-0.500,0.459} for HG. (152) 

In the case of d 23 = 7 r /2 or 37 r/ 2 , we find 

cos <5 = {0.418,0.779} for TBM; (153) 

cos 6 = {0.498, 0.761} for GRA; (154) 

cos <5 = {0.346, 0.837} for GRB; (155) 

cos 5 = {0.394, 0.906} for HG. (156) 


7.2 The Scheme with Ui3(9f 3 , Sl 3 )Ui 2 (9l 2 , ^ 12 ) (Case D2) 


We consider the following parametrisation of the PMNS matrix (see Appendix [Cj first case 
in Table fl4l): 


u = U 13 ( 613 . <5 1 e 3)Ri2(0i2) J Pi(^R23(02 O 3)^i3(^ 3 )Co , Pi( 6 ) = diag(e^, 1,1). (157) 

A sum rule for cos 5 is obtained in the cases of either 0 ? 23 = Aur, k = 0,1,2, or @ 33 = qn/2, 
q = 0,1, 2, 3,4. For the general form of U we find for the absolute value of the element U^: 


U„ 2 


cos 6 \ 2 cos 9 


231 ) 


(158) 


which in each of the two limits indicated above is fixed because | cos $ 12 ! can be expressed 
in terms of the PMNS neutrino mixing angles. This can be seen from the following relation, 
which is obtained using the expressions for 11 /^ 31 2 in the standard parametrisation of the 
PMNS matrix U and in the parametrisation given in eq. (157): 


cos 2 (9 13 sin 2 0-23 = \ ~ e lS sin 612 sin + cos 0 12 cos sin d 23 | 2 


(159) 


Equating the expression for | t/p, 2 1 given in eq. (158) with the one in the standard parametri¬ 
sation, we find 


cos 5 = 


cos 2 @23 cos 2 @12 + sin 2 @12 sin 2 @13 sin 2 @23 — cos 2 @12 cos 2 @23 
sin 2@23 sin @12 cos @12 sin @13 


(160) 


7.3 The Scheme with Ci 2(@} 2 5 ^ 12 )^ 23(^235 ^ 23 ) (Case D3) 

We consider the following parametrisation of the PMNS matrix (see AppendixjCj second case 
in Table 14): 


U = U 12 ( 6 { 2 , 6 e 12 )R 23 (d 23 )P 2 ( 6 )R l 3 (d° 13 )R 12 (d° 12 )Q 0 , P 2 (5) = diag(l, e i 6 , 1 ) 


(161) 


A sum rule for cos <5 can be derived in the cases of either 0 33 = kir, k = 0,1,2, or 6\ 2 — 
qir/2 , q = 0,1, 2, 3,4. Indeed, the relation cos 2 @13 cos 2 @23 = cos 2 @23 cos 2 (which can be 
obtained from the expressions for the eleme nt U T 3 of the PMNS matrix U in the standard 
parametrisation and in the one given in eq. (161)), allows us to express cos 2 @23 hr terms of 


the known product cos^ @13 cos^ 6 23 and the parameter cos^ @ 33 which, in principle, is fixed by 
the symmetries Gf and G u . We have also 


I U T 2 1 = | e 1 ^ cos 0° l2 sin @23 + cos @23 sin 0\ 2 sin 1 . 


(162 
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In the limits of either 0 ° 3 = kn, k = 0,1,2, or 0 \ 2 = <?vr/2, q = 0,1, 2, 3,4, \U T 2 \ does not 
depend on 5 and is also fixed. This makes it possible to derive a sum rule for cos A In the 
general case, cos 6 is a function of 5: 

COS <5 = ------ 3 —— COS 2 0?n (COS 2 0?o — COS 2 013 COS 2 023) 

sin 2012 sin 2023 sin 0i3 cos 2 0^ L V ’ 

— cos 2 0 i 2 sin 2 023 cos 2 0 33 + cos 2 023 (cos 2 0 i 3 sin 2 0 32 sin 2 0 33 — sin 2 0 i 2 sin 2 0 i 3 cos 2 0 i 3 ) 


+ K COS 6 COS 013 COS 023 Sin 2032 si 11 013 (cos 2 013 - cos 2 013 cos 2 0 23 ) 


(163) 


where k = 1 if 623 belongs to the first or third quadrant, and k = —1 otherwise. For 0 ° 3 = 0 


the sum rule reduces to the one derived in 1111 and discussed in detail in 11,13,14 


7.4 The Scheme with Ui3(9l 3 , 5l 3 )U23(0 23 , S 23 ) (Case D4) 

We consider the following parametrisation of the PMNS matrix (see AppendixjCj second case 
in Table 14): 


u = U 13 ( 0 e 13 , 5 t 3 )R 2 3 ( 023 )P 2 (S)Ri 3 ( 0 °i 3 )Ri 2 (O°i 2 )Qo , P 2 (S) = diag(l, e a , 1). (164) 

In this case a sum rule for cos <5 exists provided either 0^3 = kit, k = 0,1, 2, or 0° 2 = qn/ 2, 
q = 0,1, 2, 3,4. This follows from the relation \Ufj , 3 \ 2 = cos 2 0i3 sin 2 0 23 = cos 2 0i 3 sin 2 023 and 
the expression for | |: 

\U^ 2 \ = | e 1 ^ cos 0i 2 cos 023 — sin 023 sin0° 2 sin0i31. (165) 

The sum rule of interest for cos 5 reads 

2 


cos 5 = — 


cos 2 012 (cos 2 013 — cos 2 013 sin 2 023) 


sin 2012 sin 2023 sin @13 cos 2 0^ 3 
— cos 2 0 i 2 cos 2 023 cos 2 0 ^ + sin 2 023 (cos 2 0 i 3 sin 2 0 i 2 sin 2 0° 3 — sin 2 0 i 2 sin 2 0 i 3 cos 2 0 i 3 ) 


— k cos 5 cos 0 i 3 sin 023 sin 20 ° 2 sin 0 ° 3 (cos 2 0i 3 — cos 2 0 i 3 sin 2 023 ) 


(166) 


where k = 1 if 0 23 belongs to the first or third quadrant, and k = — 1 other wise. As in the 
previous case, cos 5 is a function of 5. For 0^3 = 0 the sum rule in eq. (166) reduces to the 
one derived in 1141. 


7.5 The Scheme with ^ 23 (^ 23 ; ^ 23 )^ 13 (^ 13 ) ^ 13 ) (Case D5) 

In this case we consider the following parametrisation of the PMNS matrix (see Appendix [Cj 
third case in Table 14): 


u = U23(0 2318 h)Ri 30 is)Pi(t)R 23 (Oh)Ri 2 (O°i 2 )Qo , Pi(S) = diag(e^ 1 , 1 ) 
We find that: 

sin 2 0 i 3 = \U e3 \ 2 = cos 2 0 23 sin 2 0 i3 , 


Sin (723 = . 


\ U ^\ _ • 2 


sin 012 = 


l-|C e3 | 2 

\Ue2\ 2 

l-|C e3 | 2 


— sin 023($13) $ 23 ) ^ 23 ’ $ 23 ) 
= sin 2 012(013; 5, $ 12 ) $ 23 ) • 


(167) 

(168) 

(169) 

(170) 
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Since, as can be shown, | U ll2 | is a function of the parameters $23> ^23’ ^13> ^12 and @ 23 , and 

6 * 13 , and cos 8 can be extracted from eqs. (168) and (170), respectively, it might be possible to 
find a sum rule for cos <5 in the case of fixed < 5 2 3 - Since in this case the analytical expression 
of cos 5 in terms of 5 23 is rather complicated, we do not present it here. Note that imposing 
either 0\ 2 = 0 or @ 2 3 = 0 is not enough to fix the value of cos <5. Even in the case of fixed 
S 23 it follows from eqs. (169) and (170) that for any given value of <5 23 , cos 5 can take several 
values. This can be understood, e.g., from eq. (170) which allows to fix cos <5, but not sin 5. 
This ambiguity, in particular, leads to multiple solutions for cos 6. In Fig. [ 6 ] we show these 
solutions in the cases of the TBM, GRA, GRB and HG symmetry forms. We remind that 
for these forms @ 23 = — 7 t /4 and @ 32 = sin~ 1 (l/v / 3) (TBM), @ 32 = sin -1 (1/^2 + r ) (GRA), 
r = (l+y / 5)/2 being the golden ratio, 0\ 2 = sin - 1 (V3 — r/2) (GRB), and 0 32 = 7 r /6 (HG). We 
assume @13 to lie in the firs t qu adrant. The solid lines correspond to 5 = cos _1 (cos <5), where 
cos 5 is the solution of eq. (170), while the dashed lines correspond to 5 = 2n — cos -1 (cos 6). 
Multiple lines reflect the fact that eq. (169) for @ 23 has several solutions. We note that Fig. [ 6 ] 
does not change in the case of @13 belonging to the third quadrant, while for @ 1 3 lying in the 
second or fourth quadrant the solid and dashed lines interchange. For <5 23 = 0 or n, we find 


cos 6 = {-0.114,0.114} for TBM; (171) 

cos 5 = {-0.289,0.289} for GRA; (172) 

cos 5 = {-0.200,0.200} for GRB; (173) 

cos 6 = {-0.476,0.476} for HG. (174) 


It is worth noting that in the scheme under consideration the values of 5 23 in a vicinity of 
tt/2 ( 37 t/ 2 ) do not provide physical values of cos 5 (see Fig. [ 6 ]). 


7.6 The Scheme with ^12(^12)^12)^13(^13)^13) (Case D 6 ) 

It is convenient to consider the following parametrisation of the PMNS matrix U (see Ap¬ 
pendix [Cj third case in Table 14): 

u = u 12 ( 0 f 2 , < 5 1 e 2 )Ri 3 ( 0 i 3 )Hi(< 5 )R 2 3 ( 02 ° 3 )^i 2 (^ 2 )Qo , Pi(S) = diag(e^, 1 , 1 ). (175) 


We find that a sum rule for cos 5 can be derived if either 0\ 2 = W2, q = 0,1, 2, 3,4, or 
@ 23 = kir, k = 0,1, 2. Indeed, the relation | t/ r 3 1 2 = cos 2 @i 3 cos 2 0 23 = cos 2 @i 3 cos 2 @ 23 , allows 
us to determine cos 2 @i 3 in terms of the known quantity cos 2 @i 3 cos 2 @ 23 and the parameter 
cos 2 @ 23 , which is fixed once Gf and G u are fixed. Further, we have 


\U t2 \ = |sin @ 32 sin@i 3 + cos @i 3 cos @ 32 sin@ 23 1 , (176) 

where the only unconstrained parameter is the phase 6. In the cases indicated above with 
either @^ 2 = qn/ 2, q = 0,1, 2, 3,4, or @ 23 = kir, k = 0,1, 2, the absolute value of the element 
U t2 does not depend on 5, which in turn allows a sum rule for cos <5 to be derived. In general, 
cos 5 is a function of 5: 


cos 6 = 


sin 2 @12 (cos 2 @23 — cos 2 @13 cos 2 @ 23 ) 


sin 2 @i 2 sin 2@ 23 sin @13 cos 2 @ 23 
— cos 2 @12 sin 2 @ 23 cos 2 @ 23 + cos 2 @ 23 (cos 2 @13 cos 2 @i 2 sin 2 @ 23 — sin 2 @i 2 sin 2 @13 cos 2 @ 23 ) 


+ k cos 5 cos @13 cos @23 sin 2 @i 2 sin @23 (cos 2 @23 — COS 2 @13 COS 2 @ 23 ) 


(177) 
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Figure 6: Dependence of cos S on <5f 3 in the cases of the TBM, GRA, GRB and HG symmetry 
forms. The mixing parameters sin 2 $ 12 , sin 2 623 and sin 2 do have been fixed to their best fit 
values for the NO neutrino mass spectrum quoted in eqs. (Ji]) - The angle di 3 is assumed 
to belong to the fir st qu adrant. The solid lines correspond to 5 = cos -1 (cos <5), where cos 8 is 
the solution of eq. (170), while the dashed lines correspond to 5 = — cos _1 (cos <5). See text 

for further details. 


where k = 1 if $13 belongs to the first or third quadrant, and k = — 1 otherwise. In this 
the sum rule for cost) has been derived first in 14 assuming d 33 = 0, but as we can see 
result holds also for any fixed value of @ 33 , since the parametrisation given in eq. (175) 
the corresponding one in |14| are the same after a redefinition of the parameters. 

The sum rules derived in Section [7] are summarised in Table 0 


case 

this 

and 
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Case 

Parametrisation of U 

Sum rule for cos 6 


D2 

C/l3(@f3^f 3 )«12(012) J Pl(5)^23(02 O 3)^13(^3)Qo 

cos 2 @23 cos 2 @12 + sin 2 @12 sin 2 @13 sin 2 @23 — cos 2 @12 cos 2 @ 23 

sin 2@23 sin @12 cos @12 sin @13 

D3 

u 12 (dt 2 ,5f 2 )R23(e 23 )p 2 (d)R 13 (ei 3 )R 12 (ei 2 )Q 0 

2 

COS 2 @12 (cos 2 @13 — cos 2 @13 COS 2 @23) 

sin 2 @i 2 sin 2@23 sin @13 cos 2 @i 3 


— cos 2 @12 sin 2 @23 cos 2 @13 + cos 2 @23 (cos 2 @13 sin 2 @12 sin 2 @i 3 — sin 2 @12 sin 2 @13 cos 2 @i 3 ) 
+re cos S cos @13 cos @23 sin 2 @i 2 sin @13 (cos 2 @13 — COS 2 @13 COS 2 @23) 2 


D4 


Ui3(ot 3 ,6f 3 )R 23 (e 23 )p 2 (8)Ri 3 (e o 13 )Ri2(e°i 2 )Q 0 


om zju \ ^ oin Zo w ^ 

— cos 2 @12 cos 2 @23 cos 2 @13 + sin 2 @23 (cos 2 @13 sin 2 @12 sin 2 @i 3 


0111 vi‘z * 


— K COS 


<5 cos @13 sin @23 sin 2 @i 2 sin @i 3 (cos 2 @i 3 — cos 2 @13 sin 2 @23) 


D6 


C/l2(@?2,<5f2)i?13(013)^l(5) J R23(^3)^12(^ 2 )Qo ■ 

sin z(7i2 sin z(723 

— cos 2 @12 sin 2 @23 cos 2 @ 23 + cos 


5 aiiit/13 


i 2 @23 (cos 2 @13 cos 2 @12 sin 2 @23 — sin 2 @12 sin 2 @13 cos 2 @ 23 ) 
+k cos 5 cos @13 cos @23 sin 2 0\ 2 sin ^23 (cos 2 @ 23 — cos 2 @13 cos 2 @23) : 


Table 7 : Summary of the sum rules for cosi5 in the case of fully broken G e under the assumption that the matrix U e consists of two 
complex rotation matrices. The parameter re = 1 if the corresponding hat angle belongs to the first or third quadrant, and re = — 1 
otherwise. The cases D3 and D4 have been analysed for @i 3 = 0 in 111.14]. In the case D6 the sum rule for cos S has been derived first 


in 14] assuming @i 3 = 0, but this result holds also for any fixed value of @i 3 . See text for further details. 





















8 The Case of Fully Broken G u 


When the discrete flavour symmetry Gf is fully broken in the neutrino sector, the matrix U v 
is unconstrained and includes, in general, three complex rotations and three phases, i.e., three 
angle and six CPV phase parameters. It is impossible to derive predictions for the mixing 
angles and CPV phases in the PMNS matrix in this case. Therefore we will consider in this 
section forms of U u corresponding to one of the rotation angle parameters being equal to zero. 
Some of these forms of U v correspond to a class of models of neutrino mass generation or 
phenomenological studies (see, e.g., [48]) and lead, in particular, to sum rules for cos 5. Since 
in this case Gf is fully broken in the neutrino sector, the Z 2 x Z 2 symmetry of the Majorana 
mass term does arise accidentally. Therefore the matrix U v is not constrained by the symmetry 
group Gf. We give in Table 14 in Appendix [C| the most general paranretrisations of U under 
the assumption that for fully broken G v one rotation angle vanishes in the matrix U v . 


8.1 The Scheme with ^12)^13(^13) ^13) (Case El) 

It proves convenient to consider the following parametrisation of the PMNS matrix U in this 
case (see Appendix [Cj fourth case in Table 14): 


U = R 23 (e° 23 )R 13 ( 0 l 3 )P 1 ( 6 )R 1 2 ( 0 12 )U 13 (e l { 3 , S 33 )Qo , Pi(tf) = diag(e^, 1 , 1 ). (178) 


Consider first the case of 9 ^ 3 = 0. In this case the phase <5 is unphysical. Comparing this 
parametrisation of U with the standard parametrisation, we find: 


sin 2 0 i 3 = | t 7 e 31 2 = sin 2 0\ 3 cos 2 #12 


sin 2 023 = 


IM 5 


l-\U e3 \ 2 COS 2 013 


[ sin 2 023 cos 2 0 i 3 + cos 2 023 sin 2 0 33 sin 2 9\2 


1 


— — sin 205 q sin 20 ^ sin 0 i 2 cos 6 


■ 2 n W 

sm 0 i 2 = 


J 13j > 


From the ratio 


_ sin 2 0 i 2 

l-|C e3 | 2 COS 2 013 ' 


U r2 


u. 


/x2 


— tan 2 0 2 3 , 


(179) 

(180) 
(181) 

(182) 


we get the following sum rule for cos 5: 


cos 5 = -— [cos 2003 sin 2 0 i 3 + (sin 2 023 — sin 2 0? 3 ) (cot 2 0 i 2 — sin 2 013 )] . (183) 

sm 2023 sm 0 i 3 

Substituting the best fit values of the neutrino mixing angles for the NO neutrino mass 
spectrum and the value of 9% 3 = — 7 r/ 4 , which corresponds to the TBM, BM, GRA, GRB and 
HG symmetry forms, we obtain cos <5 = 0.616. We note that in the considered scheme the 
predictions for cos <5 are all the same for the symmetry forms mentioned above, since these 
forms are characterised by different values of the angle 0\ 2 i which has been absorbed by the 
free parameter 0 i 2 - This “degeneracy” can be lifted in specific models in which the value of 
0^2 is fixed. Using the best fit values and the requirement | cosdj < 1, we find that the allowed 
values of sin 2 9 23 belong to the following interval: 0.338 < sin 2 9 23 < 0.538. 
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In order to give the general result for cos <5 in the case of 7 ^ 0, we use the expression 
for sin 2 0\2 for non-zero 0 J 3 : 


sin 2 6 »i 2 = 


I u t 


e21 


cos 2 0 f 3 sin 2 0 


12 


1 - 11^31 2 COS 2 O 13 

Employing this relation in the expression for |t/ r2 | 2 , we get 

2 


(184) 


cos 5 = — 


cos 2 $23 (cos 2 — sin 2 9 12 cos 2 # 13 ) 


sin 2012 sin 2023 sin 0i 3 cos 2 0 13 L 

— cos 2 0 i 2 cos 2 023 cos 2 0\ 3 + sin 2 0 i 2 (cos 2 0 i 3 sin 2 0\ 3 sin 2 02 3 — sin 2 0 i 3 sin 2 023 cos 2 0 ° 3 ) 

— k cos 5 sin 0 i 2 cos 0 i 3 sin 0i 3 sin 26 23 (cos 2 0 ° 3 — sin 2 0 i 2 cos 2 013 ) 


(185) 


where k = 1 if 0 i 2 belongs to the first or third quadrant, and k = — 1 otherwise. 


Similar to the cases C 2 , C5, C7 and C9 analysed in subsections |5.2| , [54)1 |5 . 7| and |5 . 9[ cos S 
is a function of the known neutrino mixing angles 0 i 2 , 0 i 3 and 023 , of the angles 9f 3 and 0 ^ 
fixed by Gf and the assumed symmetry breaking pattern, as well as of the phase parameter 5 
of the scheme. Predictions for cos 5 can be obtained if 5 is fixed by additional considerations 
of, e.g., GCP invariance, symmetries, etc. 

For 0 ^ 3 = kir, k = 0,1, 2, and/or 0 23 = k'n/2, k' = 0,1, 2, 3,4, cos 5 does not depend on 5 


and k. In the first case the expression in eq. (185) reduces to the sum rule given in eq. (183). 


8.2 The Scheme with ^12(^12; ^12)^23(^23) ^23) (Case E 2 ) 

In this case it is convenient to use another possible parametrisation of the PMNS matrix, the 


fourth case in Table 14 given in Appendix [Cj Namely, 


u = R 23 (e° 23 )R 13 (e o 13 )P 1 (d)R l 20 1 2 )U 2 3 ^ 23 , S 23 )Qo , Pi(S) = di a g(e iS , 1,1) 

Consider first the possibility of 0i 3 = 0. Under this assumption we find: 

sin 2 0i 3 = \U e3 \ 2 = sin 2 9 23 sin 2 0 i 2 , 

1 


sin 023 = 


l ^ 3 | 2 


l-|U e3 | 2 COS 2 013 


[ sin 2 0 : 


23 


cos 


7 23 


+ cos 02 3 sin 0 23 cos 0 


23 ' 


7 12 


1 


+ - sin 29 23 sin 29 23 cos 0 i 2 cos S 23 \ , 


sin 2 0 i 2 = 


|U e2 | 2 


cos 2 02 3 sin 2 0 i 2 


l-|U e3 | 2 COS 2 013 

The sum rule of interest for cos 5 can be derived in this case using the ratio 

2 


U< 


r 1 


u, 


m 1 


= tan 2 0 ° 


23 


(186) 

(187) 

(188) 

(189) 

(190) 


We get 
cos 6 = 


cot 012 


sin 2023 sin 0 i 3 


[cos 202 3 sin 2 0 i 3 + (sin 2 0 2 3 — sin 2 0 23 ) (tan 2 0i 2 — sin 2 @ 13 )] . (191) 
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This sum rule can be formally obtained from the r.h.s. of eq. (183) by interchanging tan $12 


and cot 0 i 2 and by multiplying it by (—1). Substituting the best fit values of the neutrino 
mixing angles for the NO neutrino mass spectrum and the value of $33 = —tt/ 4, we get 
cos 6 = —0.262. Using the best fit values and the requirement | cos<5| < 1, we find that the 
allowed values of sin 2 $23 belong to the following interval: 0.227 < sin 2 $23 < 0.659. 

In order to find a general result for cos <5 for arbitrary fixed 0 7 ^ 0, we use the following 
relation: 

cos 2 012 COS 2 013 = cos 2 012 cos 2 0i 3 , (192) 

which follows from the expressions for | C/ e i | 2 in the standard parametrisation and in the 


parametrisation given in eq. (186). With the help of this relation, using |J7^i|, we get 


cos 5 = 


cos 2 023 (cos 2 013 — cos 2 012 COS 2 013) 


sin 2012 sin 2023 sin 0 i 3 cos 2 0° 3 
— sin 2 0 i 2 cos 2 023 cos 2 0 i 3 + cos 2 0 i 2 (cos 2 0 i 3 sin 2 0 i 3 sin 2 0 ^ — sin 2 0 i 3 sin 2 023 cos 2 0 i 3 ) 


+ re cos S cos 0 i 2 cos 0 i 3 sin 013 sin 2023 (cos 2 013 - cos 2 012 COS 2 01 3 ) 


(193) 


where re = 1 if 0 i 2 belongs to the first or third quadrant, and re = — 1 otherwise. Also in this 
case cos <5 is a function of the unconstrained phase parameter 5 of the scheme. Predictions 
for cos 5 can be obtained if S is fixed by additional considerations (e.g., GCP invariance, 
symmetries, etc.). 

As like in the case El, for 0 ^ 3 = kn, k = 0,1, 2, and/or 0% 3 = k'n /2, k ' = 0,1, 2, 3,4, cos 5 
does not depend on 5 and re. For 0J 3 = 0 , 7 r, 27r, the sum rule in eq. (193) coincides with the 


sum rule given in eq. (191) 


8.3 The Scheme with ^23(^235 ^23)^12(^12; ^12) (Case E 3 ) 


The convenient parametrisation for U to use in this case is that of the fifth case in Table 14 
given in Appendix [C) 

C = i?13(01 O 3)^12(^2)^2^)i?23(02 3 )Ci2(0r 2 ,^2)Qo, P 2 (S) = diag(l, e iS , 1 ) . 

We find that: 

sin 2 0 i 3 = \U e3 \ 2 = sin 2 013(023,5, 0 i 2 , 0i 3 ), (194) 

, ,2 _2/io 2 n 

sin 2 023 = 1 M ' 3 


sin 2 0 i 2 = 


|£V 2 

i-|c e3 | 2 

We2? 

i-|c e3 | 2 


COS 2 013 

= sin 2 012(023, 6,0^,5^, 6° 12 ,9 ° 13 ). 


(195) 

(196) 


However, a sum rule for cos 6 cannot be obtained because cos 5 turns out to depend, in 
particular, on 6 ” 2 which is an unconstrained phase parameter of the scheme considered, which 
can be seen from the expression for | U tJi 1 1: 


| U^\ | = | cos 0 i 2 sin 0 i 2 + e l ^ +s " 2 ' > cos 023 cos 0 J 2 sin 0/ 2 1 • 


(197) 


The situation here is analogous to the cases analysed in subs ectio ns 7.1 and 7.5. Namely, 
considering a certain residual symmetry group G e , from eq. (195) we find that sin 2 023 is 
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fixed. Then, cos <5 is fixed (up to a sign) by eq. (194). Hence, 0 1 2 can be expressed in terms 


of 51 2 by virtue of eq. (196). Thus, numerical predictions for cos <5 can be obtained if 5 \ 2 is 
fixed. 

8.4 The Scheme with U 23 {6 23 , ^ 23 )^ 13 (^ 13 ; ^ 13 ) (Case E4) 

Employing the parametrisation for U given in Appendix [Cj namely the fifth case in Table 14 

U = Ris(ei 3 )Ru(ei 2 )P 2 (5)R 23 (§ 23 )U l3 (ei 3 , S v 13 )Q q , P 2 (6) = diag(l, e rf , 1), 

we find that cos <5 is a function of d 23 , 9 and the PMNS mixing angles. Therefore, cos 5 can 
be determined only in those cases when 0 23 is fixed. Using the result 

\U e2 \ 2 1 


sin 0i2 = 


l-|u e3 | 2 COS 2 013 


cos^ 023 cos- 0 i 3 sin 2 0 i 2 + sin 2 023 sin 2 0i 3 


1 


— - cos 5 sin 2023 sin 20^ sin 0i 2 


(198) 


we find these cases to be, for example: i) 0^ 2 = 0 , tv, leading to the relation sin 2 0 i 2 cos 2 0 i 3 = 
sin 2 023 sin 2 0 J 3 , ii) 9\ 3 = 0, tv, implying sin 2 0 i 2 cos 2 0 i 3 = cos 2 023 sin 2 0 J 2 , hi) 0 f 3 = tv/2 , 
37t/ 2, giving sin 2 0 i 2 cos 2 0 i 3 = sin 2 023. For this reason we give cos 5 as a function of the 
angle 0 23 . Namely, the sum rule of interest, which is obtained using \U^, 2 \ = | cos 023 cos0i 2 |, 
reads 

, cos 2 012 cos 2 023 + sin 2 012 sin 2 013 sin 2 023 - cos 2 023 COS 2 0?2 /-> nn \ 

cos 5 =---------— . (199) 

sin 2023 sin 0 i 2 cos 0 i 2 sin 0 i 3 

The dependence of cos 5 on Gf is realised via the values of the angles 0\ 2 anf i # 13 - 
8.5 The Scheme with Ui 3 (Q 1 3 , ^13)^12(^12; ^12) (Case E 5 ) 

The parametrisation for the PMNS matrix U employed by us in this subsection is the sixth 
case in Table [14] given in Appendix [Cj 

U = R23(0° 23 )R 12 (e o 12 )P 1 (5)R 13 (9 13 )U 12 (91 2 ,51 2 )Q 0 , Pi(5) = diag(e iS , 1,1). 

We find that: 

sin 2 0 i 3 = |U e3 | 2 = cos 2 0 i 2 sin 2 0 i3 , 


sin- 023 = 


sin- 0 i 2 = 


I U, 


li3| 


l-|C e3 | 2 


— sill 2 023(013, $12) ^23) > 

= sin 2 0 12 (0 13 ,5,O1 2 ,51 2 ,0° 12 ). 


( 200 ) 

( 201 ) 

( 202 ) 


i-|c e3 | 2 

However, a sum rule for cos 5 cannot be obtained because cos 5 turns out to depend, in 
particular, on 5\ 2 which is an unconstrained phase parameter of the scheme considered. This 
can be seen, e.g., from the expression for | U tl 1 1: 

\Uf_til = | cos 0i2(e^ sin0i2 cos 023 cos 0i3 + sin0i3 sin 02 3 ) + e^ 12 cos 0i2 cos 023 sin0i21• (203) 


Similarly to t he case analysed in subsection 8.3, for a certain residual symmetry grou p G e , 
from eq. (200) we find that sin 2 0i3 is fixed. Then, cos 5 is fixed (up to a sign) by eq. (201), and 


so the angle 0^ can be expressed in terms of 5\ 2 by virtue of eq. (202). Therefore, numerical 
predictions for cos 5 can be obtained if b\ 2 is fixed. 
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8.6 The Scheme with Ui3(0/ 3 , 5± 3 )U 23 (9% 3 , ^ 23 ) (Case E6) 


The parametrisation of the PMNS matrix U utilised by us in the present subsection is that 

given in Appendix [Cj 


of the sixth case in Table 14 


u = R23(0° 23 )R 12 (d° 12 )P 1 (d)R 13 (e 13 )U 23 (e^ 5% 3 )Qo , P^S) = diag(e rf , 1 ,1). 


A sum rule and predictions for cos <5 can be derived in the cases of either 0 23 = qn/ 2, q = 
0,1, 2, 3,4, or 61 2 = fc7r, k = 0,1, 2. Indeed, using the relation 

| Cel | 2 = COS 2 012 cos 2 013 = COS 2 013 COS 2 0i 2 , (204) 

we can express cos 2 0i3 in terms of the product of PMNS neutrino mixing parameters cos 2 0i2 
cos 2 0i3 and, the fixed by Gf parameter, cos 2 0^2. The sum rule of interest for cos <5 can be 
derived, e.g., from the expression for the absolute value of the element U^\: 

| tv I = | e - ^ cos 0i3 cos 023 sin 0^2 + sin 0i3 sin 0231, (205) 


since in any of the two limits indicated above, 0 23 = qi r/2, q = 0,1, 2, 3,4, or 0 32 = kir, 
k = 0,1, 2, \U^i\ does not depend on 6 . In fact, it is given only in terms of the known PMNS 
neutrino mixing parameters and an angle (either 0 23 or 0° o ) whi ch is fixed by the symmetry 
Gf. In the general case, cos 5 is a function of <5. Using eqs. (204) and (205[), we get 


005,5 = 5 in 20 1 . 2 sin 20 2 3 Sin 0 13 C o 5 ^; 2 

— sin 2 0 i 2 cos 2 023 cos 2 0 i 2 + cos 2 0 i 2 (cos 2 0 i 3 sin 2 0^2 cos 2 023 — sin 2 0 i 3 sin 2 d 23 cos 
+ re cos 5 cos 0i2 cos 0i3 sin 012 sin 2023 (cos 2 012 - cos 2 012 cos 2 01 3 ) : 


2 0°u) 

(206) 


where re = 1 if 0i3 lies in the first or third quadrant, and re = — 1 otherwise. For 0° 2 = kn, 
k = 0,1,2, and/or 0 23 = k'n/2, k! = 0,1,2, 3,4, cos 6 does not depend on 6 and re. 

The sum rules derived in Section [8] are summarised in Table [H 
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Case Parametrisation of U Sum rule for cos 5 

El ^ 23 ( 023 )^ 13 (^ 3)^1 (£)iM012)tfl3(0i3,$i 3 )Qo • 9 , . -n 0.0 COS 2 0^3 (cos 2 0^3 Sill 2 012 COS 2 01 

Sill ZC 712 sin ZC 723 Sin (713 cos (733 

— cos 2 0i2 cos 2 023 cos 2 0° 3 + sin 2 0i2 (cos 2 0i3 sin 2 0^3 sin 2 0? 2:i — sin 2 

— k cos 5 sin 0i2 cos 0i3 sin 0 ^ 3 sin 20 £ 3 (cos 2 — sin 2 0i2 cos 2 013) 5 

3) 

013 sill 2 023 COS 2 033) 

E 2 R2m 3 )Rim 3 )Pimn(6n)UoMf>^Z 3 )Qo M . M . , ,„ 0 

sin Zu\2 sin z(723 sin (713 cos (733 

— sin 2 012 COS 2 023 COS 2 033 + cos 

+K COS S COS 012 COS 013 sin 033 sil 

COS 2 023 (cos 2 033 — COS 2 012 COS 2 013 

s 2 0i2 (cos 2 0i3 sin 2 033 sin 2 d 23 — sin 2 

11 2023 (cos 2 033 — COS 2 012 COS 2 @13) 2 

) 

013 sin 2 023 cos 2 03 3 ) 

El R,i(0° )P,(S)R»(ff»)U,i( 0 1 ' snon cos20 12 COS 2 023 +Sin 2 012 Sin 2 013 Sin 2 023 -COS 2 023 COS 2 0^2 


E 4 Ifi 3 ( 0 i 3 )^i 2 ( 0 i 2 )^(djlf 23 ( 023 jci 3 ( 0 i 3 ^i 3 )yo sin 20 23 sin 0 12 cos 0 i2 sin 013 

sin .1(712 sin 1(723 sm (713 cos (732 

— sin 2 012 COS 2 023 COS 2 032 + cos 

+K COS S COS 012 COS 013 sin 032 si] 

sin 2 023 (cos 2 012 — COS 2 012 cos 2 013 ) 

i 2 0i2 (cos 2 0i3 sin 2 032 cos 2 023 — sin 2 

11 2023 (cos 2 032 — COS 2 012 COS 2 @13) 2 

1 

013 sin 2 023 COS 2 012 ) 


Table 8 : Summary of the sum rules for cos <5 in the case of fully broken G v under the assumption that the matrix [/„ consists of two 
complex rotation matrices. The parameter k = 1 if the corresponding hat angle belongs to the first or third quadrant, and k = — 1 
otherwise. See text for further details. 




9 Summary of the Predictions for Gf = A 4 (T'), S 4 and 


In this section we summarise the numerical results obtained in the cases of the discrete flavour 


diagonalisation of the corresponding group elements which lead to physical values of cos 5 
and phenomenologically viable results for the “standard” mixing angles 6 12 , $13 and $ 23 - In 
the cases when the standard mixing angles are not fixed by the schemes in Tables 0 HU we 
use their best fit values for the NO spectrum quoted in eqs. © @. For the cases in the 

tables marked with an asterisk, physical values of cos <5, i.e., | cos<5| < 1 , cannot be obtained 
employing the best fit values of the neutrino mixing angles 6 12 , $13 and # 23 , but they can be 
achieved for values of the relevant mixing parameters allowed at 3a. Note that unphysical 
values of cos <5, | cos <5| > 1 , occur when the relations between the parameters of the scheme 
and the standard parametrisation mixing angles cannot be fulfilled for given values of sin 2 $ 12 , 
sin 2 0i3 and sin 2 $ 23 - Indeed the parameter space of sin 2 $ 12 , sin 2 $13 and sin 2 @23 is reduced 
by these constraints coming from the schemes. 


symmetry groups A 4 (T 7 ), S 4 and A 5 , which have been already discussed in subsections 3.4 


4.4 and 5.10 In Tables IP we give the values of the fixed angles, obtained from the 


(G e , G u ) = ( Z 3 , Z 2 ) cos 6 sin 2 0 12 


B1 (sin 2 0l 2 , sin 2 0% 3 ) = (1/3,1/2) 0.570 0.341 


Table 9: The phenomenologically viable case for the symmetry group A 4 . The values of cos 6 
and sin 2 612 predicted by the scheme Bl, which refers to the corresponding parametrisation in 
Tables [ 3 ] and [HJ have been obtained using the best fit values of sin 2 #13 and sin 2 d 23 for the NO 
spectrum quoted in eqs. © ©. 

For the symmetry group A 4 we find that the residual symmetries 

• ( G e , G u ) = ( Z 2 , Z 2 ) in the cases Cl - C9; 

• ( G e , G v ) = (Z 3 , Z 2 ) in the cases B2 and B3; 

• (G e , G u ) = ( Z 2 x Z 2 , Z 2 ) in the cases Bl, B2 and B3; 

• (G e , G v ) = ( Z 2 , Z 3 ) or ( Z 2 , Z 2 x Z 2 ) in the cases Al, A2 and A3 

do not provide phenomenologically viable results for cos 5 and/or the standard mixing angles. 
It is worth noticing that the predicted value of sin 2 di 2 = 0.341 in Table 0 is within the 
2 a allowed range. Varying sin 2 $ 13 , which enters into the expression for sin 2 di 2 , within its 
respective 3a allowed range for the NO neutrino mass spectrum, we find 0.339 < sin 2 8\ 2 < 
0.343. 

For the symmetry group S 4 we find that the residual symmetries 

• ( G e , Gy) = (Z 2 , Z 2 ) in the cases C 6 and C9; 

• ( G e , Gy) = (Z 3 , Z 2 ) in the case B3; 

• (G e , Gy) = (Z 4 , Z 2 ) or (Z 2 x Z 2 , Z 2 ) in the cases B2 and B3; 

• ( G e , Gy) = ( Z 2 , Z 3 ) in the cases Al, A2 and A3; 
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(G e , G v ) = (Z2, Z4) or (Z2, Z2 X Z2) in the case A3 


do not provide phenomenologically viable results for cosh and/or for the standard mixing 
angles. 

The cases in Table 10 marked with an asterisk are discussed below. Firstly, using the best 
fit values of sin 2 0\2 and sin 2 #13 we get a physical value of cos h in the case C3 for the minimal 
value of sin 2 6*23 = 0.562, for which cosh = —0.996. For C 8 with sin 2 9% 3 = 1/2 and 3/4, using 
the best fit values of the neutrino mixing angles for the NO spectrum, we have cosh = —1.53 
and 2.04, respectively. The physical values of cosh can be obtained, using, e.g., the values of 
sin 2 #23 = 0.380 and 0.543, for which cosh = —0.995 and 0.997, respectively. In the parts of 
the 3cr allowed range of sin 2 $ 23 , 0.374 < sin 2 023 < 0.380 and 0.543 < sin 2 023 < 0.641, we 
have —0.938 > cosh > —0.995 and 0.997 > cosh > 0.045, respectively. Secondly, in the case 
B1 we obtain cos h = —0.990 employing the best fit value of sin 2 0 i 3 and the maximal value of 
sin 2 023 = 0.419. Finally, utilising the best fit value of sin 2 0 i 3 , we get physical values of cosh 
in the cases A1 and A2 for the minimal value of sin 2 0 i 2 = 0.348, for which cosh = —0.993 
and 0.993, respectively. Note that for the cases in which sin 2 023 is fixed, the predicted values 
are within the corresponding 2 cr range, while in the cases in which sin 2 0 i 2 is fixed, the values 
of sin 2 0 i 2 = 0 .341 and 0.31 1 are within 2(7 and 1 ( 7 , respectively. The value of sm~ 0 i 2 — 0.256 
lies slightly outside the current 3er allowed range. 

For the symmetry group A 5 we find that the residual symmetries 


(Ge, 

Gy) 

= (Z 2 ,Z 2 ) 

in the 

cases 

C 2 , C 6 and C7; 

(Ge, 

Gy) 

= (Z 3 , Z 2 ) 

in the 

cases 

B2 and B3; 

(Ge, 

Gy) 

= (Z 5 ,Z 2 ) 

in the 

case 

B3; 


(Ge, 

Gy) 

= (Z2 x Z2, Z2) in the 

cases 

B1 and B3; 

(Ge, 

Gy) 

= (Z 2 ,Z 3 ) 

or ( Z 2 

,z 5 ) 

in the 

case A3; 

(Ge, 

Gy) 

= (Z 2 , Z2 X Z 2 ) in the 

cases 

Al, A 2 and A3 


do not provide phenomenologically viable results for cosh and/or for the standard mixing 
angles 0 12 , 013 and 0 2 3 - 

We will describe next the cases in Table [TT] marked with an asterisk, apart from those which 
have also been found for Gf = £4 and discussed earlier. Using the best fit values of sin 2 0 i 2 and 
sin 2 0 i 3 we get a physical value of cos h in the case C4 for the minimal value of sin 2 623 = 0.487, 
for which cosh = —0.997. Instead using the best fit values of sin 2 0 i 3 and sin 2 623 one gets the 
physical values of cosh = — 1 for the maximal value of sin 2 0 i 2 = 0.277. Employing the best 
fit value of sin 2 0 i 3 we find a physical value of cos h in the case B2 with residual symmetries 
(G e , Gy) = (Z2 x Z2, Z2) for the minimal value of sin 2 623 = 0.518, for which cosh = —0.996. 
Similarly for the cases A 1 and A 2 with residual symmetries (G e ,G u ) = (Z 2 ,Z§), the values 
of cosh = —0.992 and 0.992 are obtained using the minimal value of sin 2 0 i 2 = 0.321. 

The values of sin 2 9°- in Table lTj used to compute cos h and sin 2 0,y are the following 
ones: l/(4r 2 ) = 0.0955, (3 - r)/4 = 0.3455, 1/(2 + r) = 0.2764, 1/(4 + 2 r) = 0.1382, 
1/(3 + 2 r) = 0.1604, 1/(3 + 3r) ^ 0.1273, 2/(4r 2 - r) ^ 0.2259, r/( 6 r - 6 ) ^ 0.4363, 
( 6 r - 4)/(10r - 3 ) = 0.4331, (1 - r )/(8 - 6 r) = 0.3618. 
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(G e ,G v ) = (Z 2 ,Z 2 ) 

cos 6 

sin 2 Oij 

Cl sin 2 #23 = 1/4 

-0.806 

not fixed 

C2 sin 2 0% 3 = 1/2 

not fixed 

sin 2 #23 = 0.512 

C3 sin 2 91 3 = 1/4 

- 1 * 

not fixed 

C4 sin 2 # 3 2 = 1/4 

0.992 

not fixed 

C5 sin 2 #32 = 1/4 

not fixed 

sin 2 #12 = 0.256 

C7 sin 2 #23 = 1/2 

not fixed 

sin 2 #23 = 0.488 

C 8 sin 2 #£3 = { 1 / 2 , 3/4} 

{—1*, 1*} 

not fixed 

(G e , G v ) = (Z 3 ,Z 2 ) 

cos 5 

sin 2 #12 

B1 (sin 2 # 32 , sin 2 # 23 ) = (1/3, 1 / 2 ) 

0.570 

0.341 

B2 (sin 2 # 32 , sin 2 # 33 ) = (1/6,1/5) 

-0.269 

0.317 

(G e ,G u ) = (Z 4 ,Z 2 ), (Z 2 x Z 2 ,Z 2 ) 

cos 5 

sin 2 #12 

B1 (sin 2 # 32 , sin 2 # 23 ) = (1/4,1/3) 

- 1 * 

0.256 

(G e , G u ) = (Z 2 , Z 4 ), (Z 2 ,Z 2 xZ 2 ) 

cos 6 

sin 2 #23 

Al (sin 2 # 33 , sin 2 # 23 ) = (1/3,1/4) 

-1* 

0.488 

A2 (sin 2 # 3 2 , sin 2 # 23 ) = (1/2,1/2) 

1 * 

0.512 


Table 10: The phenomenologically viable cases for the symmetry group S4. The values 
of cos S and sin 2 # 12 or sin 2 # 23 predicted by the schemes Al, A 2 , etc., which refer to the 
corresponding parametrisations in Tables i © have been obtained using the best fit values 
for the NO spectrum of the other two (not fixed) neutrino mixing parameters (sin 2 0 13 and 
sin 2 623, or sin 2 O12 and sin 2 $13) quoted in eqs. ^ |8]). In the cases marked with an asterisk, 

physical values of cos S cannot be obtained employing the best fit values of the mixing angles, 
but are possible for values of the relevant neutrino mixing parameters lying in their respective 
3 cr allowed intervals. See text for further details. 
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(' G e ,G u ) = (Z 2 ,Z 2 ) 

cos 6 

sin 2 Oij 

Cl sin 2 $23 = 1/4 

-0.806 

not fixed 

C3 sin 2 = 0.0955, 1/4 

0 . 688 , - 1 * 

not fixed 

C4 sin 2 0 ° 2 = 0.0955, 1/4 

-1*, 0.992 

not fixed 

C5 sin 2 0\ 2 = 1/4 

not fixed 

sin 2 612 = 0.256 

C8 sin 2 $23 = 3/4 

1 * 

not fixed 

C9 sin 2 0° 12 = 0.3455 

not fixed 

sin 2 0\2 = 0.330 

(G e ,G„) = (Z 3 ,Z 2 ) 

cos 5 

sin 2 0i2 

B1 (sin 2 0° 12 , sin 2 6? a ) = (1/3, 1 / 2 ) 

0.570 

0.341 

(G e , G v ) = (Z 5 , Z 2 ) 

cos 5 

sin 2 0i2 

B1 (sin 2 0J 2 , sin 2 0? 2:i ) = (0.2764,1/2) 

0.655 

0.283 

B2 (sin 2 0° 2 , sin 2 0° 3 ) = (0.1382,0.1604) 

-0.229 

0.259 

(G e ,G u ) = (Z 2 xZ 2 ,Z 2 ) 

cos 6 

sin 2 0 i 2 

B2 (sin 2 0° l2 , sin 2 0f 3 ) = (0.0955,0.2764) 

- 1 * 

0.330 

(1/4,0.1273) 

0.805 

0.330 

(G e , G v ) = (Z 2 ,Z 3 ) 

cos 6 

sin 2 023 

Al (sin 2 0J 3 , sin 2 6? 23 ) = (0.2259,0.4363) 

0.716 

0.553 

A2 (sin 2 0 32 , sin 2 0 23 ) = (0.2259,0.4363) 

-0.716 

0.447 

(G e ,G v ) = (Z 2 ,Z 5 ) 

cos 6 

sin 2 023 

Al (sin 2 0 33 , sin 2 0 23 ) = (0.4331,0.3618) 

-1* 

0.630 

A2 (sin 2 ^ 2 , sin 2 ^ 3 ) = (0.4331,0.3618) 

l* 

0.370 


Table 11: The phenomenologically viable cases for the symmetry group A5. The values 
of cos S and sin 2 0 i 2 or sin 2 $23 predicted by the schemes Al, A 2 , etc., which refer to the 
corresponding parametrisations in Tables i i have been obtained using the best fit values 
of the other standard mixing angles for the NO spectrum quoted in eqs. {6} In the 

cases marked with an asterisk, the predicted values of cos 6 , obtained for the best fit values 
of the neutrino mixing angles $12, 0 13 and $23, are unphysical; physical values of cos <5 can be 
obtained for values of the neutrino mixing parameters sin 2 $12, sin 2 $13 and sin 2 $23 lying in 
their respective 3 a allowed intervals. See text for further details. 
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10 Conclusions 


In the present article we have employed the discrete symmetry approach to understanding 
the observed pattern of 3 -neutrino mixing and, within this approach, have derived sum rules 
and predictions for the Dirac phase 6 present in the PMNS neutrino mixing matrix U . The 
approach is based on the assumption of the existence at some energy scale of a (lepton) flavour 
symmetry corresponding to a non-Abelian discrete group Gf. The flavour symmetry group 
Gf can be broken, in general, to different “residual symmetry” subgroups G e and G u of the 
charged lepton and neutrino mass terms, respectively. Given Gf , typically there are more 
than one (but still a finite number of) possible residual symmetries G e and G u . The residual 
symmetries can constrain the forms of the 3 x 3 unitary matrices U e and U v , which diagonalise 
the charged lepton and neutrino mass matrices, and the product of which represents the PMNS 
neutrino mixing matrix U, U = ul U v . Thus, by constraining the form of the matrices U e 
and U u , the residual symmetries constrain also the form of the PMNS matrix U . This can 
lead, in particular, to a correlation between the values of the PMNS neutrino mixing angles 
#12, #i3 and @23) which have been determined experimentally with a rather good precision, 
and the value of the cosine of the Dirac CP violation phase <5 present in U, i.e., to a “sum 
rule” for cos A The sum rule for cos <5 thus obtained depends on residual symmetries G e and 
G v and in some cases can involve, in addition to @12, @13 and @23) parameters which cannot 
be constrained even when Gf is fixed. For a given fixed Gf, unambiguous predictions for the 
value of cos 6 can be derived in the cases when, apart from the parameters determined by Gf 
(and G e and G u ), only @12, @13 and @23 enter into the expression for the respective sum rule. 

In the present article we have derived sum rules for cos 5 considering the following discrete 
residual symmetries. 1) 17 g — ^2 and — Zyi, u 7 * 2 or Z^i x Z^^, ?r, an ^ 2 (Section [3]); 
ii) G e = Z n , n > 2 or Z n x Z. m , n,m> 2 and G v = Z2 (Section]!]); iii) G e = Zo and G u = Z% 
(Section [5]); iv) G e is fully broken and G u = Z n , n > 2 or Z n x Z m , n,m > 2 (Section [ 7 ]); 
and v) G e = Z n , n > 2 or Z n x Z m , n,m > 2 and G u is fully broken (Section [ 8 ]). The sum 
rules are summarised in Tables ii 0 and [8] For given G e and G u , the sum rules for cos <5 
we have derived are exact, within the approach employed, and are valid, in particular, for 
any Gf containing G e and G v as subgroups. We have identified the cases when the value of 
cos 5 cannot be determined, or cannot be uniquely determined, from the sum rule without 
making additional assumptions on unconstrained parameters (cases A 3 in Section [ 3 ] and B 3 
in Section i (see also Table [ 3 ]); cases C 2 , C 5 , C 6 , C 7 and C 9 in Section [5] (see also Table [I]); 
the cases discussed in Sections [ 7 ] and [ 8 ]) . In the majority of the phenomenologically viable 
cases we have considered the value of cos 5 can be unambiguously predicted once the flavour 
symmetry Gf is fixed. In certain cases of fixed Gf, G e and G v , correlations between the 
values of some of the measured neutrino mixing parameters sin 2 @12, sin 2 @13 and sin 2 @23, 
are predicted, and/or the values of some of these parameters, typically of sin 2 @12 or sin 2 @23, 
are fixed. These correlations and “predictions” are summarised in Tables [ 5 ] and [ 6 ] We have 
found that a relatively large number of these cases are not phenomenologically viable, i.e., 
they lead to results which are not compatible with the existing data on neutrino mixing. 
We have derived predictions for cos <5 for the flavour symmetry groups Gf = S4, A4, T' 
and A5 using the best fit values of sin 2 @12, sin 2 @13 and sin 2 @23, when cos 5 is unambiguously 
determined by the corresponding sum rule. We have presented the predictions for cos <5 only in 
the phenomenologically viable cases, i.e., when the measured values of the 3 -neutrino mixing 
parameters sin 2 @12, sin 2 @13 and sin 2 @23, taking into account their respective 3 cr uncertainties, 
are successfully reproduced. These predictions, together with the predictions for the value 
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of one of the mixing parameters sin 2 @12 and sin 2 023 ; in the cases when it is fixed by the 
symmetries, are summarised in Tables 0-0 

The results derived in the present study show, in particular, that with the accumulation 
of more precise data on the PMNS neutrino mixing parameters sin 2 0]2, sin 2 613 and sin 2 023; 
and with the measurement of the Dirac phase 5 present in the neutrino mixing matrix U, it 
will be possible to critically test the predictions of the current phenomenologically viable the¬ 
ories, models and schemes of neutrino mixing based on different non-Abelian discrete (lepton) 
flavour symmetries Gf and sets of their non-trivial subgroups of residual symmetries G e and 
G u , operative respectively in the charged lepton and neutrino sectors, and thus critically test 
the discrete symmetry approach to understanding the observed pattern of neutrino mixing. 
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A Appendix: The Discrete Groups A 4 , T', S 4 and A* 


A 4 is the symmetry group of even permutations of four objects (see, e.g., !)• it is isomorphic 
to the tetrahedral symmetry group, i.e., the group of rotational symmetries of a regular 
tetrahedron. As such it can be defined in terms of two generators S and T, satisfying S 2 = 


T 3 = (ST ) 3 = 1. In this work, we choose to work in the Altarelli-Feruglio basis 145 for the 
3-dinrensional representation of the S and T generators, see Table [12] 


The group T' is the double covering group of A 4 (see, e.g., |2|j), which can be defined in 
terms of two generators S and T through the algebraic relations: R 2 = T 3 = (ST ) 3 = 1, 
RT = TR, where R = S 2 . We use the basis for the 3-dimensional representation of the 


generators S and T from 30 , summarised in Table 12 Since we restrict ourselves to the triplet 


representation for the LH charged lepton and neutrino fields, there is no way to distinguish 
T' from A 4 30 15 Note that matrices representing S and T in Table 12 for A 4 , are related 
with those for T' by the following redefinition S —> TST 2 , T —>■ T 2 , where S and T before 


(after) the arrows are the matrices presented in Table 12 for A 4 (T' 


S 4 is the group of permutations of four objects, i.e., the rotational symmetry group of a 
cube (see, e.g., |2|). It can be defined in terms of three generators S, T and U, satisfying 


46 


S 2 = T 3 = U 2 = (ST ) 3 = ( SU ) 2 = ( TU ) 2 = ( STJJ ) 4 = 1. We employ for the 3-dimensional 


representation of the S, T and U generators the basis given in |46 and summarised in Table 12 
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As it was also shown in |46|, this basis is equivalent to the basis widely used in the literature 


A$ is the group of even permutations of five objects (see, e.g., [2]), i.e., the rotational 
symmetry group of an icosahedron, which can be defined in terms of two generators S and 


5 It is worth noting that A4 is not a subgroup of T'. 
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Group 3-dimensional representation of the generators 


A 4 


T 


S 4 


A 5 



T = 


T = 


0 0 \ 
u; 2 0 

0 u) 

0 0 \ 
w 0 
0 u 2 


S 


S 



0 

1 

0 



T 


1 

7E 



-V2 

—r 

l/r 



T 


/ i —y/2i —A 

\V2 0 vM 

\ i \/2i —i) 



0 

-1 

0 


0 

0 


(l 0 

0 p 

\0 0 



Table 12: 3-dimensional representation of the generators of A 4 , T', S 4 and A 5 . We have 
defined u> = e 2 ™/ 3 , r = (1 + -\/5)/2 and p = 


T, satisfying S 2 = T° = (ST ) 3 = 1 . We employ the basis defined in 47 , which for the 


3-dimensional representation of the generators S and T is summarised in Table 12 


We conclude this appendix by noting that a list of the Abelian subgroups of A4, T', 
and A 5 can be found in 


£4 


49 , 17|, 46 and 47 , respectively. 


B Appendix: Parametrisations of a 3 x 3 Unitary Matrix 

Parametrisations of a 3x3 unitary matrix W (see, e.g., |50[|52|) can be obtained, e.g., from one 
of the six permutations of a product of three complex rotations and diagonal phase matrices, 
e.g., as follows: 

W = ^ 2 ^ 3 W = *4*2*3 Uij U kl U rs , (207) 

where we have assumed ij ^ kl ^ rs. It is worth noticing that sometimes it is convenient to 
use the parametrisations of W of the following form: 


II — Uij U k i Uij . (208) 

As shown in |50|, the number of distinctive parametrisations of a CKM-like matrix is nine. 


We have defined the phase matrices T, in eq. (16) and the complex rotation matrix in the i-j 
plane Uij = Uij(9ij,5ij) in eq. ( [T7| ). The latter can be always parametrised as a product of 
diagonal phase matrices and the rotation matrix Rij = Rij(9ij) = Uij(9ij, 0 ), i.e., 

Uij = Pi(8)* Rij Pi(6) = Pj(—5)* R^ Pji-6 ), (209) 

where P t (S) are diagonal matrices defined as follows: 

Pi(S) = diag(e i( 5 , 1 , 1 ), P 2 (d) = diag(l,e i( 5 ,l), P 3 (S) = diag(l, 1 , e iS ) . ( 210 ) 
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( 211 ) 


Defining Pij(a,/3 ) as a product Pij(a,/3 ) = Pi(a)Pj(f3), the following relation holds: 


Uij (%, 8ij) Pij (a, (5) = Pij (a, /3)Uij (%, 5',-), 


with 5'ij = Sij + a — /3. 


Starting from the general parametrisation of W in eq. ( 207[ ) and the relation in eq. (211), 
we find convenient parametrisations for W. They are summarised in Table 13. The parametri- 
sations of the matrix 6/°(0f 2 , ^ 13 , $231 {<%}) defined in Section [ 2 ] have been obtained from 
Table 13 after a redefinition of the phases {5^}. For example, in the first case when 
U°(e° 12 ,e° 13 ,e° 23 ,{S° kl }) is represented by the product U 12 (0° 2 , 8 l 2 ) u 23 ( 0 23 ,5^ 3 )U 13 {91 3 ,5l 3 ) the 


Case 

Initial form of W 

Final parametrisation of W 

Al 

U 12 U 23 U 13 

^ 12(^131 ^ 23 )^ 12 (^ 12 , <^12 — di 3 + 823 ) R 23 P 13 P 12 (5i3, ^ 2 , 3 ) 

A2 

U 13 U 23 U 12 

F**3($12, —^ 23 )^ 13 ( 6 * 13 , 5i3 — £12 — S 23 ) P-23 7612 Pi 3 (51 2 , —5 2 3) 

A3 

U 23 U 13 U 12 

F23(^12, ^ 13 ) 6 / 23 ( 6 * 23 , <5^23 + < 5 l 2 — 5 i 3 )Ri 3 Ri 2 P 2 3 { 5 l 2 , ^ 13 ) 

B1 

$ 

CO 

Cl 

to 

00 

-Pl 3 (£i 2 , —^23)-R23-Rl2t/l3(013, <^13 — 8\2 — S23)Pl3{$12, - 823 ) 

B2 

U 13 U 12 U 23 

P 23 (5l2, ^ 13 )-6^13-6^126/23 (^23, <^23 + 5l2 — <^ 13 ) 

B3 

U 23 U 13 U 12 

^12(^13, ^23)^23-^136/12(^12, <5i2 — 5l3 + ^23)-Pl2(^13, ^23) 

Cl 

U 12 U 23 U 13 

P 3 (^ 23 ) 6 / 12 (^ 12 , 5i2)-6?23//l3(6*13, <^13 ~ 8 23 )P 3 (8 2 3) 

C2 

U 13 U 23 U 12 

P 3 (<^23)//l3(013, ^13 — <>23)-P23f/l2(#12, ^ 12)^3 (<523) 

C3 

U 12 U 13 U 23 

-F 3 (^ 13 ) 6 / 12 (^ 12 , 8i2)Rl3U23(0 2 3, &23 ~ 813 ) (^ 13 ) 

C4 

U 13 U 12 U 23 

P 2 (^12)//l3(013, <^13)-Pl2t/23(^23, <^23 + 8 l 2 )P 2 (^ 12 ) 

C5 

U 23 U 12 U 13 

-f* 2 (^ 12 ) 1/23(^23, 823 + 5l2)76l 2 6/13 ( 013 , <5l3)P|(<5l2) 

C 6 

U 23 U 13 U 12 

-P3(<5l3)6/23(023, 823 ~ <^13)-Pl36/l2(012, 8 \ 2 )P 3 (<^ 13 ) 

C7 

U 12 U 23 U 12 

(^23)6/l2(012, <^12)76236/l2(012, <512)7^(523) 

C 8 

U 13 U 23 U 13 

P 2 (^23)6/13(013, 5l3)76236/i3(0i3, 5i3)P2(523) 

C9 

U 23 U 12 U 23 

P] (5l2)6/23(023, 523)76i26/23(023, 523)7^1(512) 


Table 13: Equivalent parametrisations of W obtained using the result in eq. (211), which 
allows us to find the convenient form of the matrix U°(0^ 2l 9 ° 3 , 9 23 , {5^}) defined in Section[2j 


following redefinition is used: 5 32 ~ 5J 3 + 8 23 ~^ ^12 • 

The product of two complex rotations in the i-j plane can always be written as 

= Pij(P, —a)R ij (O ij )Pi(a - (5) = Pj(-a - P)R ij (9 ij )P ij (a, (3) (212) 
= Pij(a, -p)R ij {9 ij )P j {P -a) = Pi(a + P)R ij {9 ij )P ij {~P, -a ), 
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where we have defined the angle $j~ as 


sin On = I s%c\,e + cg.-sf.-e 


„6 -»<& 


'y 


y y 


y y 


and the phases a, f3 as 


a 


= arg [ 


c f. c b _ 


] , P = arg [s$C(j e + 44e **] , 


h i<5t. i 


(213) 

(214) 


with = sin<9“- &) and c“- b) = cos $“- 6) . 


C Appendix: The Case of Fully Broken G e or G v 

In the case when the group G e is fully broken and G u = Z n , n > 2 or Z n x Z m , n,m > 2, 
there are cases in which one can express cos <5 as a function of $\ 2l $ 13 , $23 and $f 2) $13 > $ 23 - 
In the cases 

i ) u\ = £43(13) ($23(13), $J3(13) )££l2($l2> ^ 12 )) 
ii) *4 = ££l2(13) ($12(13) ’ ^12(13))^23 ($23 > ^2 3 ) ) 
iii) ££] = ££ 23 ( 12 ) ($ 23 ( 12 ) > ^ 23 ( 12 )) ££13 ($13 > ^i 3 )j 
we choose for convenience, respectively: 

i) ££°($ 12 ; $13j $23) ^ 12 ) = ££ 12 ($ 12 ) ^12)££23($23)-££l3($13) : 

ii) C£°($ 12 ) $13) $23) ^ 23 ) = ££23($23) ^23)^13($ 13 )^ 12 ($ 12 ) ) 

iii) ££°($ 12 ; $ 13 > $23) ^ 13 ) = ££l3 ($( 3 ) ° 3 )-£*23 ($ 23 )£*"12 ($( 2 ) • 

The possib le pa ra metri sations of U presented in Table [l4| can be obtained from i), ii) and iii) 
using eqs. (212) - (214). The angles $f), $ ; y and the phases Sf^, 5 are free parameters. It can 
be seen from Table H that if one of the fixed angles turns out to be zero, the number of free 


parameters reduces from four to three. The same situation happens if one of the two free 
phases is fixed. Thus, in some of these cases a sum rule for cos <5 can be derived. 

In the case when the group G e = Z n , n > 2 or Z n x Z m , n,m> 2 and G v is fully broken, 
we consider the following forms of the matrix U v , 

iv) U u = ££ 12 ($ 12 , ^12)££l3(23)($13(23)) 5 13(23))<9o, 

v) U v = ££23 ($23) ^23) ££l2(13) ($12(13)) ^12(13) )Q0) 
vi) U u = ££l 3 ($r 3 ) ^3)££l2(23)($r2(23)) ^2(23))Qo, 

for which we choose for convenience, respectively: 


iv) ££° ($12) $13: $23) ^12 ) 

V) ££°($i 2 , $13; $23: ^23 ) 
Vi) ££°($12) $13) $23) ^13) 


= R2 3 (e° 23 )Ri3(e° 13 )u 12 (e° 12 ,6° 12 ), 
= Ri3(0° 13 )R 12 (9 o 12 )U 23 (9 o 23 ,S o 23 ), 
= R23(0° 23 )R 12 (9° 12 )U 13 (9° 13 ,6° 13 ). 


The para metris ati ons o f U in the cases iv), v) and vi) presented in Table 14 have been obtained 


with eqs. (212) - (214). The angles $f-, $,, and the phases 5’/-. 5 are free parameters. It can 
be seen from Table |14 that if one of the fixed angles turns out to be zero, the number of free 
parameters reduces from four to three. The same situation happens if one of the two free 
phases is fixed. Thus, in some of these cases a sum rule for cos <5 can be derived. 
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^(0f2,0f3.0l3.{^i}) 

Parametrisation of U for fully broken G e 

^23(13) (023(13)’ ^23(13))^ 12 (012) ^ 12 ) 
^12(13)^2(13)’ <5l2(13)) [/ 23(02 e 3^2 e 3) 
^23(12)(0 2 e 3(12)^2 e 3 ( 12))^3(0?3’^3) 

^23(13)^3(13)’ 5 2 e 3(13))^12(012)A(<5)i?23(02 O 3)^13(0f3)Qo 

G2(13)^2(13)’ 5 12(13))^3(023)P2(5) J Rl3(^ 3 )^12(0?2)Qo 

^23(12)(02 e 3(12)’ 5 2 e 3(12))^13(013)A(^ J R2 3 (02°3)^12(01 O 2)Qo 

^(0?2’013’023’{%}) 

Parametrisation of U for fully broken G u 

U\ 2 {9\ 2 i ^12)G3(23)(013(23)’ ^13(23)) 

^23(02^3’ ^3)^12(13)(0^2(13).'^2(13)) 
^13 (013 > ^13)^12(23)(012(23)’ ^12(23)) 

i?23(02 O 3)^13(0?3)^l(^^12(012)Ui3(23)(0r 3 (23)’^3(23))Qo 

i?13(0l O 3)i?12(^2)^2(^i?23(023)G2(13)(0r2(13)’^2(13))^O 

^23 (0° 23 ).R 12 (9° 12 )P 1 (5) R 13 (013 ) Ui 2(23 ) (012(23) ’ ^2(23) )Qo 


Table 14: Upper (lower) part. Parametrisations of U in the case of fully broken G e (G „) and 
= Z n , n > 2 or Z n x Z rn , n,m> 2 (G e = Z n , n > 2 or x Z m , n, m > 2) when U e ((/„) 
has particular forms. 


D Appendix: Results for Gf = A§ and Generalised CP 


Models with A ,5 and GCP symmetry have been recently developed by several authors |8 10 


We show that our results for the symmetry group A 5 under the same assumptions of [101 and 
the same breaking patterns reduce to the one derived in 10 . The results in eqs. (10), (11), 


(12) and (14) in 101 lead to the following phenomenologically viable cases: 


i )U = diag(l, i, -i)R 23 (0° 23 )R 12 (012) diag(l, -*, i)Ri 3 ffi 3 ) , for G e = Z 3 , G v = Z 2 , 
ii) U = diag(l, i, -i)R2z{dl 3 )R 12 {0l 2 ) diag(l, -i, i)R 13 {9i 3 ) , for G e = Z 5 , G v = Z 2 , 
hi) U = diag(l, 1, -l)R 23 {9° 23 )R l 2 (9° l2 ) diag(l, 1, -l)^^), for G e = Z 5 , G u = Z 2 , 
iv) U = Ri 3 (9° 13 )R 12 (9° 12 )R 23 (9° 23 ) diag(l, 1, -1)R 23 (0" 23 ) , for G e = Z 2 x Z 2 , G u = Z 2 , 


where we have in i) 0 \ 2 = sin 1 (l/y / 3) and #23 = — vr/4, ii) 9 \ 2 = sin 1 (l/\/2 + r) and 
= -k /4:, ii i) 91 2 = sin~ 1 (l/V2 + r) a nd 9° 23 = —tt/ 4, iv) 9° l2 = sin- 1 (l/(2r)), 9° l3 = 
sin _1 (l/-v/2 + r) and 9 23 = sin” 1 (r/\/2 + r). 

Using (sin 2 0 32 , sin 2 0f 3 , sin 2 0 23 ) = (1/3,0,1/2) in the case i), the results in eqs. (56)-(58), 
after defining 9 13 = 0 33 = 9 and setting d '13 = (5 33 = 7 r/ 2 , reduce to 

2 1 2 1 

sin 9 23 = — and cos (5 = 0. 


sin“ 0 i 3 = - sin^ 9 , sin^ 0 i 2 =-= 

3 3 — 2 sin 2 1 


Denoting *13 = ^3 = 0 and setting < 5 i 3 = <5 33 = 7 r /2 in the case ii), the results in eqs. (p 6 | - 
(58) reduce to 

sin 2 0 13 = 


sin 


1 + (1 — r ) 2 


sin 2 0 i 2 = 


1 


1 + r 2 cos 2 0 ’ 


sin 2 023 = - and cos (5 = 0. 
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The difference between the case iii) and the case ii) consists only in the phase <$13 which now 
is equal to tv, < 5 i 3 = = tv. Therefore while sin 2 $13 and sin 2 #12 remain unchanged, we find 


. 9 1 (sin 9 - y/1 + r 2 cos #) 2 . 

sin #23 = -——5-yv- and cos o = l. 

2 1 + cos^ # 


Finally, in the case iv) from eqs. (64) ( 66 ), defining 623 = #23 — #23 = $23 — ^ an d ^23 = 0, 

we get: 

1 + r(cos 2 9 — sin 29) 


sin 2 #13 = 
sin 2 #12 = 


l + (l-r)/(#) 


sin 2 #23 = 


2 1 + (1 — r)(cos 2 # + sin 2 #) 


3-(l-r)/(#) 


3 — (1 — r)/(#) 
and | cos <5| = 1 , 


4.2 


where /(#) = (sin 2 # — sin 2 #). Therefore the general results derived in Sections 4.1 and 
with the choices as in i), ii), iii) and iv) and the additional restriction of the parameters due 


to the presence of GCP allow one to find the formulae derived in 10 


E Appendix: General Statement 

In this appendix we prove the general statement that Z 2 symmetries preserved in the neutrino 
and charged lepton sectors can lead to phenomenologically viable predictions, only if their 
generators do not belong to the same Z 2 x Z 2 subgroup of the original flavour symmetry group. 
We compute the form of U° in a model independent way. Given a Z 2 x Z 2 symmetry with 
elements Z 2 X Z 2 = {1, gi, g 2 , # 3 } and a unitary matrix V such that V^g\V = diag(l, —1, —1), 
V^g 2 ~V = diag(—1,1, —1), V^g^V = diag(—1, —1,1), we consider first the case of G e = Z 2 = 
{1 ,gi} and G u = Z 2 = {1 ,gfl\ with i,j = 1,2,3 for all the cases Cl - C9 in Table |4j In the 
case Cl (C2) we find that the matrix U° reads 



(1 

0 

°\ 


(l 

0 

°\ 


U° = TV23 = 

0 

0 

1 for i = j , 

U° = 

0 

1 

0 for 1 j , 

(215) 


1° 

1 

«/ 


1° 

0 

V 



defined up to permutations of the 1st and 3rd (1st and 2nd) columns and the 1st and 2nd (1st 
and 3rd) rows. These permutations are not relevant because they correspond to a redefinition 
of the free parameters in the transformations f7i2(#f 2) ^ 12 ), Ui3(9^,5f 3 ) and phase matrices 
contributing to the Majorana phases or removed with a redefinition of the charged lepton 
fields. In the case C3 (C 6 ) we find that the matrix U° reads 



(° 

0 

1\ 


A 

0 

°\ 


U° = 7T 13 = 

0 

1 

* 

=2 

O 

u° = 

0 

1 

0 for i j , 

(216) 


l 1 

0 

0 J 


1° 

0 

V 



defined up to permutations of the 2nd and 3rd (1st and 2nd) columns and the 1st and 2nd 
(2nd and 3rd) rows. For the case C4 (C5) we find that the matrix U° reads 



(0 

1 

°\ 


(1 

0 

°\ 


U° = TV\2 = 

1 

0 

0 for i = j , 

U° = 

0 

1 

0 for i^j, 

(217) 


1° 

0 

V 


1° 

0 

V 
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defined up to permutations of the 2nd and 3rd (1st and 3rd) columns and the 1st and 3rd 
(2nd and 3rd) rows. The freedom in permuting the columns and rows as we described above 
does not have physical implications because it represents the freedom to perform a fixed U ( 2 ) 
transformation in the degenerate subspace of the generator of the corresponding Z 2 symmetry. 
For the other cases we find similar results. Namely, 


U° = diag(l, 1,1) for i = j and U° = 7 r 2 3 ( 13 ) for i / j for case C7, (218) 

U° = diag(l, 1,1) for i = j and U° = 7 r 23 ( 12 ) for i 7 ^ j for case C 8 , (219) 

U° = diag(l, 1,1) for i = j and U° = 7 r 13 ( 12 ) for i ^ j for case C9. (220) 


The cases in eqs. (215) - (220) do not lead to phenomenologically viable results because 
some of the elements of the resulting PMNS mixing matrix equal zero. The cases when 
a) G e = Z 2 x Z 2 = { 1 , 51 , 52 , 93 } and G u = Z 2 = {1, 55 }, b) G v = Z 2 x Z 2 = { 1 , 51 , 92 , 93 } 
and G e = Z 2 = { 1 , 9 ,;}, c) G e = Z 2 x Z 2 = { 1 , 51 , 52 , 93 } and G u = Z 2 x Z 2 = { 1 , 91 , 92 , 93 } 
are not phenomenologically viable as well. This can be seen trivially setting one or two of the 
free rotation angles, , 0 ^, to zero. 
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